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ABSTRACT 

Many topics are discussed in b physics including lifetimes, decay mechanisms, determi¬ 
nations of the CKM matrix elements VA and V u b |, facilities for b quark studies, neutral 
B meson mixing, rare b decays, hadronic decays and CP violation. 1 review techniques for 
finding physics beyond the Standard Model and describe the complimentarity of b decay 
measurements in elucidating new physics that could be found at higher energy machines. 
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1 Introduction: The Standard Model and B Decays 

Studies of b and c physics are focused on two main goals. The first is to look for new 
phenomena beyond the Standard Model. The second is to measure Standard Model 
parameters including CKM elements and decay constants. These lectures concern “B 
Phenomenology,” a topic so broad that it can include almost anything concerning b quark 
decay or production. 1 will cover an eclectic ensemble of topics that 1 find interesting and 
hope will be educational. 

1.1 Theoretical Basis 

are comprised of left-handed doublets con- 
singlets (Rosner 2001) 

Ur, d,R, Cr, Sr, tR, bR ( 1 ) 

e Ri h'Ri T Ri VeRi V VR1 V tR- (2) 

The gauge bosons, W ± , 7 and Z° couple to mixtures of the physical d, s and b states. 
This mixing is described by the Cabibbo-Kobayashi-Maskawa (CKM) matrix (see below) 
(Kobayashi 1973). 

The Lagrangian for charged current weak decays is 

A* = ~4= J^wi + h.c., 


The physical states of the “Standard Model” 
taining leptons and quarks and right handed 



( 3 ) 
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where 


and 



1 e L N 


f M 

J£ c = (z/ e , Z/ M , v T ) YVmns 

hL 

+ {ul, Cl, tL ) 7 ^VcKM 

Sl 


V J 


V b L ) 



1 Vud 

Vus 

v ub \ 

Vckm — 

Vcd 

V cs 

V cb 


V v td 

V ts 

V tb ) 


( 4 ) 


( 5 ) 


Multiplying the mass eigenstates (dr, s'r, bp) by the CKM matrix leads to the weak 
eigenstates. Vmns is the analogous matrix required for massive neutrinos (we will not 
discuss this matrix any further). There are nine complex CKM elements. These 18 num¬ 
bers can be reduced to four independent quantities by applying unitarity constraints and 
the fact that the phases of the quark wave functions are arbitrary. These four remain¬ 
ing numbers are fundamental constants of nature that need to be determined from 
experiment, like any other fundamental constant such as a or G. In the Wolfenstein ap¬ 
proximation the matrix is written in order A 3 for the real part and A 4 for the imaginary 
part as (Wolfenstein 1983) 


Vckm 


( 1 — A 2 /2 A AX 3 (p — ip)(l — A 2 /2) \ 

-A 1 - A 2 /2 - ipA 2 X A AX 2 (l + ipX 2 ) 

^ AA 3 (1 — p — irf) —AX 2 1 ) 


( 6 ) 


The constants A and A are determined from charged-current weak decays. The measured 
values are A = 0.2205±0.0018 and A=0.784±0.043. There are constraints on p and p from 
other measurements that we will discuss. Usually the matrix is viewed only up to order 
A 3 . To explain CP violation in the K° system the term of order A 4 in V cs is necessary. 


1.1.1 Determination of Gf 


Muons, being lighter than the lightest hadrons, must decay purely into leptons. The 
process is /i~ —> e~V e v M as shown in Figure [|. 

The total width for this decay process is given by 


r,= 


G 2 

J 1 m 5 x (phase space) x (radiative corrections) 
1927T d p 


( 7 ) 


Since T M • = h, measuring the muon lifetime gives a direct measure of Gp. 
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Figure 1. Diagram for muon decay. 


71° 

Figure 2. Semileptonic K~ decay diagram. 

1.1.2 Determination of |V^ S | 

A charged current decay diagram for strange quark decay is shown in Figure |j. Here the 
CKM element V us is present. The decay rate is given by a formula similar to equation (7]), 
with the muon mass replaced by the s-quark mass and an additional factor of | 14 s | 2 . 
Two complications arise since we are now measuring a decay process involving hadrons, 
K~ —» n 0 e~u rather than elementary constituents. One is that the s-quark mass is not 
well defined and the other is that we must make corrections for the probability that the 
u-spectator-quark indeed forms a ix° with the w-quark from the s-quark decay. These 
considerations will be discussed in greater detail in the semileptonic B decays section. 
Fortunately there are theoretical calculations that allow for a relatively precise measure¬ 
ment because they deal with hadron rather than quark masses and have good constraints 
on the form-factors; using the models we have A = V us = 0.2205 ± 0.0018. A is found by 
measuring V c b in semileptonic b decays and constraints on p and r) are found from other 
measurements. These will also be discussed later. 

1.1.3 Plan for Using Semileptonic B Decays to Determine V cb and V u b 

Semileptonic b decays arise from a similar diagram to Figure j|, where the b quark replaces 
the s quark. In this case the b quark can decay into either the c quark or the the u quark, 
so we can use these decays to determine V c b and V u & providing we have three ingredients 

1. B lifetimes 

2. Relevant B branching fractions 
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3. Theory or model to take care of hadronic physics. 

1.2 Lifetime Measurements 

The “ 6 -lifetime” was first measured at the 30 GeV e + e~ colliders PEP and PETRA where 
6 -quarks were produced via the diagrams shown in Figure [H| The measured the average 
lifetime of all 6 -hadron species. The distribution they found most useful was the “impact 
parameter,” which is the minimum distance of approach of a track from the primary 
production vertex. This distance is related to the lifetime (Atwood 1994). 

A more direct measurement would be to measure the actual decay distance L and 
the momentum of the 6 hadron. Then, since L = y(3ct, where t is the decay time of the 
individual particle (also called the proper time) (3 = p/E and 7 = E/mb, the distribution 
of decay times t can be derived Events will be distributed exponentially in t as e _f / T , with 
r being the lifetime. Uncertainty results from errors on L, momenta and contributions of 
backgrounds. 

Precision lifetimes of individual 6 -flavored hadrons have been measured at LEP where 
the production process is e + e _ —> Z° —> 66 and at CDF in 1.8 TeV pp collisions. Large 
samples of semileptonic decays have been used to determined the B° and B~ lifetimes. 
(Note that the CPT theorem guarantees that the lifetime of the anti-particles is the same 
as the particles.) The decay distributions for two semileptonic B decay channels are 
shown in Figure j|. The B° —> D*~£ + u channel has mostly signal with some background 
from B~ decays and other backgrounds as indicated in the figure. It takes a great deal of 
careful work to accurately estimate these background contributions. The clear exponential 
lifetime shapes can be seen in these plots. Some data has also been obtained using purely 
hadronic final states (Sharma 1994). 

A summary of the lifetimes of specific 6 -flavored hadrons is given in Figure (Groom 
2001). Note that the ratio of B + to B° lifetimes is 1.074±0.028, a 2.6cr difference from 
unity, bordering on significance. Also, the A 5 lifetime is much shorter than the B° lifetime. 
According to proponents of the Heavy Quark Expansion model, there should be at most 
a 10% difference between them (Bigi 1997). To understand lifetime differences we must 
first analyze hadronic decays. 


1.3 B Decay Mechanisms 

Figure [5| shows sample diagrams for B decays. Semileptonic decays are depicted in Fig¬ 
ure ^](a), when the virtual W~ materializes as a lepton-antineutrino pair. The name 
“semileptonic” is given, since there are both hadrons and leptons in the final state. The 
leptons arise from the virtual W~, while the hadrons come from the coupling of the spec¬ 
tator anti-quark with either the c or u quark from the 6 quark decay. Note that the B 
is massive enough that all three lepton species can be produced. The simple spectator 
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ALEPH 



t(ps) 



Kps) 


Figure 3. Proper time distributions of exclusive semi-leptonic decays from ALEPH. 



1.653±0.028 ps 
1.548 ±0.032 ps 
1.493 ±0.062 ps 
1.23±0.08 ps 


x(ps) 


Figure 4. Measured lifetimes of different B species. 


diagram for hadronic decays (Figure |5|(a)) occurs when the virtual W~ materializes as a 
quark-antiquark pair, rather than a lepton pair. The terminology simple spectator comes 
from viewing the decay of the b quark, while ignoring the presence of the spectator an¬ 
tiquark. If the colors of the quarks from the virtual W~ are the same as the initial b 
quark, then the color suppressed diagram, Figure |5|(c), can occur. While the amount of 
color suppression is not well understood, a good first order guess is that these modes are 
suppressed in amplitude by the color factor 1/3 and thus in rate by 1/9, with respect to 

















B Phenomenology 


9 




b) hadronic: color suppressed 



£, u, c 
— 1 -v, d, s 

c) annihilation 


b 

d 


W 


d) W exchange 


corn 

u 



the non-color suppressed spectator diagram. 

The annihilation diagram shown in Figure |^(c) occurs when the b quark and spectator 
anti-quark find themselves in the same space-time region and annihilate by coupling to a 
virtual W~. The probability of such a wave function overlap between the b and h-quarks 
is proportional to a numerical factor called fg. The decay amplitude is also proportional 
to the coupling V u b. The mixing and penguin diagrams will be discussed later. 

Each diagram contributes differently to the decay width of the individual species. 
Diagram (a) is expected to be dominant. There are even more diagrams expected for 
baryons. Currently there is no direct evidence for diagrams (c) and (d), although (c) is 
expected to occur, indeed it would be responsible for the purely leptonic decay B~ —> t~i 7. 

The semileptonic decay width, T s i, is defined as the decay rate in units of inverse 
seconds into a hadron (or hadrons) plus a lepton-antineutrino pair. (Decay rates can also 
be expressed in units of MeV by multiplying by ti.) T s ; is related to the semileptonic 
branching ratio B s i and the lifetime r as 

r s ; = B s i ■ r total = B s i/t . (8) 

The semileptonic width should be equal for all b species. This is true for D° and D + 
mesons, even though their lifetimes differ by more than a factor of two. Thus, it is 
differences in the hadronic widths among the different b species that drive the lifetime 
differences. 
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Let us now consider the case of B° and B~ lifetime differences. There is some indication 
that the B has a shorter lifetime, that would imply that there are more decay channels 
available. Figure |6| shows the color allowed and color suppressed decay diagrams for two- 
body decays into a ground-state charmed meson and a 7r - . The color suppressed diagram 
only exists for the B~. The relative rate 


r(B~ -»• d°tt~) 
t(b° -»• d+tt-) 


1.8 ±0.3 , 


( 9 ) 


and the same is true for all other similar two-body channels, such as D*p~. Thus we 
would expect, if most B decays are given by these diagrams that the B~ would have a 
shorter lifetime than the B , opposite of what the data suggests. An explanation is that 
this ratio reverses for higher multiplicity decays, but this is an interesting discrepancy 
that needs to be kept in mind. 



(q-d) 

(q-u) 



Figure 6. (left) Spectator diagram for B° —> D + ir and B —> D°tt . (right) Color 
suppressed spectator diagram for B~ —> D°tt~ only. 


2 Semileptonic B Decays 

2.1 Formalism of Exclusive Semileptonic B Decays 

The same type of semileptonic charged current decays used to find V us are used to find V c b 
and V u b. The basic diagram is shown in Figure ^(a). We can use either inclusive decays, 
where we look only at the lepton and ignore the hadronic system at the lower vertex, 
or exclusive decays where we focus on a particular single hadron. Theory currently can 
predict either the inclusive decay rate, or the exclusive decay rate when there is only a 
single hadron in the final state. 

Now let us briefly go through the mathematical formalism of semileptonic decays. We 
start with pseudoscalar B to pseudoscalar m transitions. The decay amplitude is given 
by (Grinstein 1986) (Gilman 1990) 

G 

A(B —> me~u) = where (10) 

v 2 


B l = u e Y (1 - y 5 ) v u , and 


( 11 ) 
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Up = {™\JUO)\B) = U(?)(P+P V + f-(v 2 )(P - PU (12) 

where q 2 is the four-momentum transfer squared between the B and the m, and P(p) are 
four-vectors of the B(m). H M is the most general form the hadronic matrix element can 
have. It is written in terms of the unknown / functions that are called “form-factors.” It 
turns out that the term multiplying the /_(g 2 ) form-factor is the mass of lepton squared. 
Thus for electrons and muons (but not r’s), the decay width is given by 


K 


dP si 
dq 2 


1 

2 Mg 


G ^ri/ + tf)lt where 

(13) 

(m 2 b + m 2 — q 2 ') — Am 2 M b 

(14) 


is the momentum of the particle m (with mass m) in the B rest frame. In principle, 
dT s i/dq 2 can be measured over all q 2 . Thus the shape of f+(q 2 ) can be determined 
experimentally. However, the normalization, /+(0) must be obtained from theory, for Vij 
to be measured. In other words, 


r SL oc |^| 2 |/ + (0)| 2 — / K 3 g(q 2 )dq 2 , 
Tb j 


(15) 


where g(q 2 ) = f + (q 2 )/f + ( 0). Measurements of scmileptonic B decays give the integral 
term, while the lifetimes are measured separately, allowing the product |V^-| 2 |/ + (0)| 2 to 
be experimentally determined. 

For pseudoscalar to vector transitions there are three independent form-factors whose 
shapes and normalizations must be determined (Richman 1995). 


2.2 Measurement of \V c b\ 

2.2.1 Heavy Quark Effective Theory and B —> D*l~v 

We can use exclusive B decays to find V c b coupled with “Heavy Quark Effective Theory” 
(HQET) (Isgur 1994). We start with a quick introduction to this theory. It is difficult to 
solve QCD at long distances, but it is possible at short distances. Asymptotic freedom, 
the fact that the strong coupling constant a s becomes weak in processes with large q 2 , 
allows perturbative calculations. Large distances are of the order 1/A qcd 1 fm, since 
A qcd is about 0.2 GeV. Short distances, on the other hand, are of the order of the quark 
Compton wavelength; A q ~ 1 jrriQ equals 0.04 fm for the b quark and 0.13 fm for the c 
quark. 

For hadrons, on the order of 1 fm, the light quarks are sensitive only to the heavy 
quark’s color electric held, not the flavor or spin direction. Thus, as rriq —» oo, hadronic 
systems which differ only in flavor or heavy quark spin have the same configuration of 
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their light degrees of freedom. The following two predictions follow immediately (the 
actual experimental values are shown below): 


m Ba - m Bd 

= m Da - rn n + 

(16) 

(90.2 ± 2.5) MeV 

(99.2 ± 0.5) MeV , and 


771 jj* Trig 

— 771 jj* 771 jj . 

(17) 

0.49 GeV 2 

0.55 GeV 2 . 



The agreement is quite good but not exceptional. Since the charmed quark is not 
that heavy, there is some heavy quark symmetry breaking. This must be accounted for 
in quantitative predictions, and can probably explain the discrepancies above. The basic 
idea is that if you replace a b quark with a c quark moving at the same velocity, there 
should only be small and calculable changes. 

In lowest order HQET there is only one form-factor function F which is a function of 
the Lorentz invariant four-velocity transfer u, where 

_ Ml + M|, - q 2 
2M r M fj* 

The point u equals one corresponds to the situation where the B decays to a D* which 
is at rest in the B frame. Here the “universal” form-factor function F(uj ) has the value, 
F( 1) = 1, in lowest order. This is the point in phase space where the b quark changes to 
a c quark with zero velocity transfer. The idea is to measure the decay rate at this point, 
since we know the value of the form-factor, namely unity, and then apply the hopefully 
small and hopefully well understood corrections. Although this analysis can be applied 
to B —» D£~u, overall decay rate is only 40% of D*£~v and the decay rate vanishes at to 
equals 1 much faster, making the measurement worse. 



2.2.2 Detection of B —» D*£v 

Since this is a semileptonic final state containing a missing neutrino, the decay cannot be 
identified or reconstructed by merely measuring the 4-vectors of the final state particles. 
One technique used in the past relies on evaluating the missing mass (MM) where 

MM 2 = ( Eb — Ed* — El) 2 — {~]t B — ~j? D * — ~p’e) 2 (19) 

= M b + M|,« + M 2 — 2 E B ■ ( Ed * + El) + 2 E b *Ei 

-Tf D * ■ + 2 ~f B ■ 77 d* + ~7e) • 


For experiments using e + e“ —> T(4S') —> BB, the B energy, E Bl is set equal to the beam 
energy, E beam , and all quantities are known except the angle between the B direction and 
the sum of the D* and lepton 3-vectors (the last term). A reasonable estimate of MM 2 is 
obtained by setting this term to zero. The signal for the D*£u final state should appear 
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at the neutrino mass, namely at MM 2 = 0. In an alternative technique MM 2 is set to 
zero and the angle between the B momentum and the sum of the D* and lepton 3-vectors 
is evaluated as 


cos (Ob-d*?) 


2 Eb(Eb* + Eg) — Mg — My,*^ 
2 I f B \ ( Pd- + Pt) 


( 20 ) 


where Mjj*g indicates the invariant mass of the _D*-lepton combination. 


A Monte-Carlo simulation of cos (@s.D*r) is given for the final state of interest and for 
the main background reaction in Figure [?]. For the correct final state only a few events 
are outside the “legal” region of ±1, while when there are extra pions in the final state 
the shape changes and many events are below —1. 



Figure 7. The cosine of the angle between the B momentum vector and the sum of D* 
and lepton momentum vectors for (left) the final state D*£u and (right) D*Xdv, where X 
refers to an additional pion. 


Recent CLEO data has been analyzed with such a technique. The data in two specific 
u bins is shown in Figure |||. The final result for all u bins is shown in Figure |j. The 
result is characterized by both a value for F(l)\V c b\ and a shape parameter p 2 . 


2.2.3 Evaluation of V ch Using B —» D*£is 

Figure 0 and Table [T] give recent experimental results on exclusive B —* D*£u decays. 
The CLEO results are not in particularly good agreement with the rest of the world 
including the BELLE results. 

To extract the value of \V c b\ we have to determine the corrections to F(l) that lower 
its value from unity. The corrections are of two types: quark mass, characterized as 
some coefficient times A qcd/itiq, and hard gluon, characterized as t]a■ The value of the 
form-factor can then be expressed as (Neubert 1996) 

F( 1) = tja (l + 0 ■ A Qou/mQ + c 2 • (A QCD /m Q ) 2 + ....) = r) A ( 1 + 6). (21) 

The zero coefficient in front of the 1 /mg term reflects the fact that the first order cor¬ 
rection in quark mass vanishes at u> equals one. This is called Luke’s Theorem (Luke 
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Figure 8. The signal and background contributions in two different u bins for the final 
state D*£u. 



Figure 9. The CLEO results for D*£u for both B° and B . The curve is a fit to a 
shape suggested in (Caprini 1998). 
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Table 1. Modern Determinations of F{l)\V c b\ using B —> D*£ v decays 


Experiment F(l)\V c b\ (xlO 3 ) p 2 


ALEPH (Buskulic 1997) 

33.0 ± 

2.1 ± 

1.6 

0.74 ±0.3 ±0.4 

BELLE (Tajirna 2001) 

35.4 ± 

1.9 ± 

1.9 

1.35 ±0.17 ±0.18 

CLEO (Heltsley 2001) 

42.2 ± 

1.3 ± 

1.8 

1.61 ±0.09 ±0.21 

DELPHI (Abreu 2001) 

34.5 ± 

1.4 ± 

2.5 

1.2 ±0.1 ±0.4 

OPAL(vrQ (Abbiendi 2000) 

37.9 ± 

1.3 ± 

2.4 

1.2 ±0.2 ±0.4 

OPAL (Abbiendi 2000) 

37.5 ± 

1.7 ± 

1.8 

1.4 ±0.2 ±0.2 


Average 37.8 ± 1.4 1.37± 0.13 



30 40 

F(l)IV cb 1x10-3 


Figure 10. The correlation between the slope parameter p 2 and F(l)\V c j,\. The contours 
are for a change in the fit \ 2 of one unit. 

1990). Recent estimates are 0.967±0.007 and —0.55±0.025 for 77^4 and S, respectively. 
The value predicted for F( 1) then is 0.91±0.05. This is the conclusion of the PDG re¬ 
view done by Artuso and Barberio (Artuso 2001). There has been much controversy 
surrounding the theoretical prediction of this number. In the future Lattice-Gauge The¬ 
ory calculations will presumably become accurate when unquenched. (Quenched calcu¬ 
lations are those performed without light quark loops.) Current lattice calculations give 
F(l) = 0.913lo;oi7=t0.016jlo;oi4iaoi6-aoi4) where the uncertainties come respectively, from 
statistics and fitting, matching lattice gauge theory to QCD, lattice spacing dependence, 
light quark mass effects and the quenching approximation. (Hashimoto 2001). 
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Using the Artuso-Barberio value for F( 1) we have 

|U c6 | = (41.5 ± 1.5 ±2.3) x 10“ 3 , 


( 22 ) 


where the first error is experimental and the second the error on the calculation of F( 1). 


2.2.4 \Vcb\ from Inclusive Semileptonic Decays 


The inclusive semileptonic branching ratio b —» X£v has been measured by both CLEO 
and LEP to reasonable accuracy. CLEO finds (10.49 ± 0.17 ± 0.43)%, while LEP has 
(10.56 ± 0.11 ± 0.18)%. These are not quite the same quantities as the CLEO number is 
an average over B~ and B° only, while the LEP number is a weighted average over all b 
hadron species produced in Z° decays. Thus using the LEP number one should use the 
average “6 quark” lifetime of 1.560±0.014 ps. 

Using the Heavy Quark Expansion, HQE, model (Bigi 1997) we can relate the total 
semileptonic decay rate at the quark level to V cb as 

\V c b\ 2 = h(Xi, A 2 , A) x T(b -> du) (23) 

= h(\i,\ 2 ,A)xB(b^c£u)/r b , 

where B{b —> c£u) is the inclusive semileptonic branching ratio minus a small b —> u£u 
component. (It is precisely the decay of a B meson to a lepton-antineutrino pair plus any 
charmed hadron.) r b is the lifetime of that particular meson or average lifetime of the 
combination of the 6-havored hadrons used in the analysis, suitably weighted. In HQE 
the semileptonic rate is described to order (A QCD/'^ib) 2 by the parameters: 


Ai = *?(B(v) 

b quark in the 


h v (iD) 2 h v B(v)), is the kinetic energy of the residual motion of the 
hadron 


A 2 = ~^(B( v) 


h 9 . /, 

1 L v 2 u V - T fits' u v 


quark spin to the gluon held. 


B(v))i is the chromo-magnetic coupling of the b 


• A = M B —m b + 2 j^, is the strong interaction coupling where M B is the spin averaged 
B meson mass, (. M B + 3M B «)/4. 


These parameters are further related as 

M B = m t + A-Li±lM (24) 

Zm b 

2m b 

These relations allow us to determine X 2 from the M B * - Mb mass splitting as 0.12 GeV 2 . 
The function h( Ai,A 2 ,A) can be calculated from the Heavy Quark Expansion (HQE) 
model. This involves both perturbative and non-perturbative pieces. 
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Although we will go through this example there is a disturbing aspect of assuming 
quark-hadron duality; the idea of duality is that if you integrate over enough exclusive 
charm bound states and enough phase space, the inclusive hadronic result will match the 
qnark level calculation. However, we do not know what size is the uncertainty associ¬ 
ated with the duality assumption. In fact, Isgnr said “I identify a source of ^.qcd/'^q 
corrections to the assumption of quark-hadron duality in the application of heavy quark 
methods to inclusive heavy quark decays. These corrections could substantially affect the 
accuracy of such methods in practical applications and in particular compromise their 
utility for the extraction of the CKM matrix element V cb (Isgur 1999). 

Let us move to the details of the calculation. In one implementation Ai and A are 
derived; the relationship between the inclusive b —> civ branching fraction (about 99% of 
b —> Xiv) is given as 


\V cb 


0.04111 


\B(b^X c iv) L55ps 


0.105 


n 


x (l ± 0.015 peri . ± 0.010 m6 ± 0.012^3 ) , 


-0.5 GeV 2 \ 
0.1 GeV 2 ) 


(25) 

(26) 


where /i 2 is the negative of Ai modulo QCD corrections and is taken as (0.5 ± 0.1) GeV 2 
(Bigi 1997). 

This leads to a value of 


| V cb \ = (40.7 ± 0.5 ± 2.4) x 1(T 3 


(27) 


from the LEP data alone, with a similar value from CLEO. The first error contains the 
statistical and systematic error from the experiments while the second error contains an 
estimate of the theory error from sources other than duality. 

In another implementation the parameters Ai and A are obtained from data. Here we 
use the HQE formula 


T,, = gHb£^ Q.369x 


192vr 3 
a, 

1 - 1.54— 

7r 


(28) 


1.65- 


A 

Mu 


a, 

1.087— 

7T 


0.95 


A 

Ml 


Ai Ao 

3.18^r + 0.02 


Ml 


Ml 


Determining Ai and A can be accomplished by measuring, for example, the “moments” 
of the hadronic mass produced in b —» civ decays. The first moment is defined as the de¬ 
viation from the D mass (Mu) as (M 2 - — M 2 D ) and the second moment as ( (m‘1 — ). 

It is also possible to use the first and second moments of the lepton energy distribution 
in these decays, or moments of the photon energy in the process b —> sy. In fact any 
two distributions can be used; in practice it will be critical to use all of them to try and 
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ascertain if any violations of quark-hadron duality appear and to check that terms of 
order (Aqcd/^&) 3 are not important. 

CLEO has used the first and second moments of the hadron mass in b —> civ decays. 
They find the Mx distributions by using missing energy and momentum in the event 
to define the v four-vector. Then detecting only the lepton and requiring it to have a 
momentum above 1.5 GeV/c, they calculate: 

M'x = (. E b -Et- E v f - {-f B l? v f (29) 

= Mf + — 2 Eb {Ei + E v ) + 2ff B ■ — r j? 1 f) (30) 

~ M b + Mf v — 2 E b E( m , 

where M( u is the invariant mass of the lepton-neutrino pair. The measured M\ distribu¬ 
tion is shown in Figure |lT| We do not see distinct peaks at the mass of the D and D* 
mesons because ignoring the last term in equation (|29|) causes poor resolution. This term 
must be ignored, however, using this technique because we do not know the direction of 
the B meson. 



Figure 11. The mass distribution for b —> c£u events from CLEO. Ah indicates an 
additional pion plus a D or D*. 


CLEO finds values of the first and second moments of (0.287±0.023±0.061) GeV 2 
and (0.712±0.056±0.176) GeV 4 , respectively. These lead to the determination of Ai, A 
and V c b shown in Figure 12. Later we will see a different determination using b —► sy 
(section [5.3.1|) . 


In summary the exclusive measurements of V c b are to be trusted while the inclusive 
determination, though consistent, has an unknown source of systematic error and should 
not be used now. 
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v cb rn b 



Figure 12. Constraints in Ai versus A from CLEO measurements of first and second 
hadronic moments in semileptonic decay. The darker band gives experimental uncertain¬ 
ties alone, while the lighter outer band also includes uncertainties from unknown 3rd order 
theoretical parameters. Bands of constant V cb and b quark mass, m b , are also shown, where 
the band widths represents theoretical uncertainties due to unknown 3rd order parameters. 

2.3 Measurement of \V u b\ 

This is a heavy to light quark transition where HQET cannot be used directly as in finding 
V c b■ Unfortunately the theoretical models that can be used to extract a value from the 
data do not currently give precise predictions. 

Three techniques have been used. The first measurement of V ub done by CLEO (Fulton 
1990) and subsequently confirmed by ARGUS (Albrecht 1990), used only leptons which 
were more energetic than those that could come from b —> c£~u decays. These “endpoint 
leptons” can occur b — » c background free at the T(45), because the B 1 s are almost at 
rest. The CLEO data are shown in Figure [13]. Since the lepton momentum for B —> D£u 
decays is cut off by phase space, this data provides incontrovertible evidence for b —> uiu 
decays. 

Unfortunately, there is only a small fraction of the b —> u£~v lepton spectrum that can 
be seen this way, leading to model dependent errors. The models used are either inclu¬ 
sive predictions, sums of exclusive channels, or both (Isgur 1995) (Bauer 1989) (Korner 
1989) (Melikhov 1996) (Altarclli 1982) (Ramirez 1990). The average among the models 
is \V ub /V cb \ = 0.079 ± 0.006, without a model dependent error. These models differ by at 
most 11%, making it tempting to assign a ±6% error. However, there is no quantitative 
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Lepton Momentum (GeV/c) 


Figure 13. Sum of inclusive electron and muon distributions from CLEO. The solid 
points are data taken on the peak of the T(4 S) while the open circles are data taken on 
the continuum 30 MeV below the resonance(suitably normalized). The dashed line is a fit 
to the continuum data and the solid line is the predicted curve from b —> civ dominated 
by B Dtv near the end of the allowed lepton spectrum. 

way of estimating the error. 

ALEPH (Barate 1999), L3 (Acciarri 1998) and DELPHI (Abreu 2000) isolate a class 
of events where the hadron system associated with the lepton is enriched in b —> u and 
thus depleted in b —> c. They define a likelihood that hadron tracks come from b decay by 
using a large number of variables including, vertex information, transverse momentum, 
not being a kaon etc.. Then they require the hadronic mass to be less than 1.6 GeV, 
which greatly reduces b —> c, since a completely reconstructed b —> c decay has a mass 
greater than that of the D (1.83 GeV). They then examine the lepton energy distribution, 
shown in Figure [L4] for DELPHI. 

The average of all three results is \V u b\ = (4.±0.20) x 10 -3 , resulting in a 
value for I4&/K.6 = 0.102±0.018, using \V c b\ = 0.0405±0.0025. I have several misgivings 
about this result. First of all the experiments have to understand the systematic errors 
very well. To understand semileptonic b and c decays and thus find their b uiv 
efficiency, they employ different models and Monte Carlo manifestations of these models. 
To find the error they take half the spread that different models give. This alone may be 
a serious underestimate. Secondly they use one model, the HQE model, to translate their 
measured rate to a value for 14,61 • This model assumes duality, and there are no successful 
experimental checks: The model fails on the A^ lifetime prediction. Furthermore, the 
quoted theoretical error, even in the context of the model, has been estimated by Neubert 
to be much larger at 10% (Neubert 2000). Others have questioned the effect of the hadron 
mass cut and estimate 10-20% errors due to this alone (Bauer 2001). 

It may be possible to use the spectrum of photons in b —> sy to reduce the theoretical 
error in the endpoint lepton method or to make judicious cuts in q 2 instead of hadronic 
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Figure 14. The lepton energy distribution in the B rest frame from DELPHI. The data 
have been enriched in b —> u events, and the mass of the recoiling hadronic system is 
required to be below 1.6 GeV. The points indicate data, the light shaded region, the fitted 
background and the dark shaded region, the fitted b —> ulv signal. 

mass to help reduce the theoretical errors. See (Wise 2001) for an erudite discussion of 
these points. 

The third method uses exclusive decays. CLEO has measured the decay rates for the 
exclusive final states irlu and plu (Alexander 1996). The data are shown in Figure |l^. 
The model of Korner and Schuler (KS) is ruled out by the measured ratio of p/n. Other 
models include those of (Isgur 1995) (Isgur 1989) (Wirbel 1985) (Bauer 1989) (Korner 
1988) (Mclikhov 1996) (Altarelli 1982) (Ramierz 1990). CLEO has presented an updated 
analysis for plv where they have used several different models to evaluate their efficiencies 
and extract V ub . These theoretical approaches include quark models, light cone sum rules 
(LCRS), and lattice QCD. The CLEO values are shown in Table ^]. 


Table 2 . Values of V u b using B —> pi v and some theoretical models 


Model 

v ub ( 

xl 0 ~ 3 ) 

ISGW2 (Isgur 1989) 

3.23 

± 

n 14+ 0 - 22 

u - i4 -0.29 

Beyer/Melikhov (Beyer 1998) 

3.32 

± 

f) 1 rr+0-21 
u - i J -0.30 

Ligeti/Wise (Ligeti 1996) 

2.92 

± 

o.i3;°i® 

LCSR (Ball 1998) 

3.45 

± 

o.15±8;|? 

UKQCD (Debbio 1998) 

3.32 

± 

0-14 ±8;§J 


The uncertainties in the quark model calculations (first three in the table) are guessed 
to be 25-50% in the rate. The Ligeti/Wise model uses charm data and SU(3) symmetry 
to reduce the model dependent errors. The other models estimate their errors at about 
30-50% in the rate, leading conservatively to a 25% error in | V ub \. Note that the models 
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Figure 15. The B candidate mass distributions and the signal bin lepton momentum 
spectra (insets) for the pion modes (top) and the sum of p and to (vector) modes (bottom). 
The points are the data after continuum and fake background subtractions; the dark- 
shaded, cross-hatched and unshaded histograms are b —> cX, b —> u£v feed-down and 
signal, respectively. For the 7 r (vector) modes, the light-shaded and hatched histograms 
are tt —> 7r (vector•—> vector) and vector —> 7r (it —> vector) cross-feed, respectively (charge 
final states can feed neutral and vice-versa). The histogram normalizations are from the 
nominal fit. The arrows indicate the momentum cuts. 


differ by 18%, but it would be incorrect to assume that this spread allows us to take a 
smaller error. It may be that the models share common assumptions, e.g. the shape of the 
form-factors. At this time it is prudent to assign a 25% model dependent error realizing 
that the errors in the models cannot be averaged. The fact that the models do not differ 
much allows us to comfortably assign a central value \V U (\ = ( 3 . 25 ± 0 . 14 %q; 29 ± 0 . 80 ) x 10” 3 , 
and a derived value \V u b/V c b\ = 0.08 ± 0.02 . 

Lattice QCD has predicted form-factors and resulting rates for the exclusive semilep- 
tonic final states it tv and ptv (Sachrajda 1999) in the quenched approximation. These 
calculations require the momentum of the final-state light meson to be small in order to 
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avoid discretization errors. This means we only obtain results at large values of the invari¬ 
ant four-momentum transfer squared, q 2 . Figure [16] shows the predictions of the B —> plv 
width as a function of q 2 . Note that the horizontal scale is highly zero suppressed. The 
region marked “phase space only” is not calculated but estimated using a phase space 
extrapolation from the last lattice point. 



Figure 16. The UKQCD lattice calculation for B —» piu shown as circles. The line is 
an estimate. 

The integral over the region for q 2 > 14 GeV 2 gives a rate of AT = 8.3 |1 2 ps -1 - 
GeV 2 . The CLEO measurement in the same interval gives (7.1 ± 2.4) x 10~ 5 ps _1 -GeV 2 , 
which yields a value for V ub = (2.9 ±0.5) x 10 -3 (Sachrajda 1999). Ultimately unquenched 
lattice calculations when coupled with more precise data will yield a much better value 
for V ub . 


3 Facilities for b Studies 

3.1 b Production Mechanisms 

Although most of what is known about b physics presently has been obtained from e + e - 
colliders operating either at the T (IS 1 ) or at LEP, interesting information is now appearing 
from the hadron collider experiments, CDF and DO, which were designed to look for 
considerably higher energy phenomena. The appeal of hadron colliders arises mainly 
from the large measured b cross-sections. At the FNAL collider, 1.8 TeV in the pp center- 
of-mass, the cross-section has been measured as ~100 ph, while it is expected to be about 
five times higher at the LHC (Artuso 1994). 

The different production mechanisms of b quarks at various accelerators leads to dis¬ 
similar methods of measurements. Figure 0 shows the production of B and B° mesons 
at the T(4S'), while Figure [18] shows the production mechanism of the different b species 
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at a higher energy e + e collider such as LEP. Figure p~9| shows the production mechanisms 
for a heavy b or c quark. 



B + B or 
B°B° 


Figure 17. B production at the T(4S'). 


b} B-,B°B s ,A b , S b> ... 

6 } 

Figure 18. b production in the continuum at e + e~ colliders. 




Figure 19. Heavy quark, Q, production at hadron colliders, (a) Second order in a s , 
while (b) is third order. 

At the Y(4S) the total b production cross-section is 1.05 nb. In hadron colliders 
the measured cross-section is about 100/rb. The third order diagrams appear about as 
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important as the second order diagrams and the overall theoretical calculation gives about 
1/2 of the measured value. 

3.2 Accelerators for b Physics 

Experiments on b decays started with CLEO and ARGUS using e + e _ colliders operating 
at the T(dS 1 ) and were quickly joined by the PEP and PETRA machines operating around 
30 GeV. Table || lists some of the machines used to study b quarks in the last century 
(Artuso 1994). 


Table 3. Machines used for b physics in the 20th century. The total number of bb pairs 
accumulated per experiment is also listed when known. 


Machine 

Beams 

Energy 

(GeV) 

o(b) 

b fraction 

C 

—2 —1 
cm s 

Total ff 

6 -pairs 

CESR 

e + e~ 

10.8 

1.05 nb 

0.25 

1.3xl0 33 

9.8 

X 

10 6 

DORIS 

e + e _ 

10.8 

1.05 nb 

0.23 

~ 10 31 

0.4 

X 

10 6 

PEP 

e + e~ 

29 

0.4 nb 

0.09 

3.2xl0 31 




PETRA 

e + e~ 

35 

0.3 nb 

0.09 

1.7xl0 31 




LEP 

e + e _ 

91 

9.2 nb 

0.22 

2.4xl0 31 

1.8 

X 

10 6 

SLC 

e + e _ 

91 

9.2 nb 

0.22 

3.0xl0 30 

8.8 

X 

10 4 

TEVATRON 

VP 

1800 

100 /ib 

0.002 

3xl0 31 





In the year 2000 the PEP II and KEK-B storage ring accelerators began operation. 
These machines have separate e~ and e + magnet rings so they can operate at asymmetric 
energies; PEP II has beam energies of 9.0 GeV and 3.1 GeV, while KEK-B has energies 
of 8.0 GeV and 3.5 GeV. This allows the B meson to move with velocity (3 ~ 0.6, which 
turns out to be very important in measurements of CP violation, since time integrated CP 
violation via mixing must be exactly zero due to the C odd nature of the T(4S'). These 
machines also make very high luminosities. Current and future machines for B physics 
are listed in Table |4j. The CDF and DO experiments will continue at the Tevatron with 
higher luminosities. CDF has already made significant contributions including studies of 
b production, lifetimes and the discovery of the B c meson (Abe 1998). BTeV and LHCb 
will go into operation around 2007 with much larger event rates. The CMS and ATLAS 
experiments at the LHC will also contribute to b physics especially in the early stages 
when the luminosity will be relatively low; at design luminosity these experiments have 
an average of 23 interactions per crossing making finding of detached vertices difficult. 
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Table 4. Machines in use or approved for dedicated b physics experiments. 


Machine 

Exp. 

Beam 

Energy 

(GeV) 

a(b) 

b 

fraction 

£ (Design) 

—2 —1 
cm s 

Interactions 

per crossing 

PEP II 

BABAR 

e + e“ 

10.8 

1.05 nb 

1/4 

3xl0 33 (f) 

< 1 

KEK-B 

BELLE 

e + e - 

10.8 

1.05 nb 

1/4 

o 

CO 

< 1 

HERA 

HERA-b 

pN 

800 

10 nb 

2•10” 6 


4 

Tevatron 

BTeV 

VP 

2000 

100 pb 

1/500 

2xl0 32 

2 

LHC 

LHCb 

VP 

14000 

500 pb 

1/160 

2xl0 32 

0.6 


f Machine has already exceeded design luminosity. 


3.3 e + e Detectors 

Most experiments at e + e _ storage rings look quite similar. CLEO II, shown in Figure pl| , 
was the first detector to have both an excellent tracking system and an excellent elec¬ 
tromagnetic calorimeter. Starting from the inside there is a thin beryllium beam pipe 



Figure 20. Electrons view of the CLEO II detector. 

surrounded by a silicon vertex detector; this detector measures positions very accurately 
on the order of 10 /im. Then there is a wire drift chamber whose main function is to mea¬ 
sure the curving trajectories of particles in the 1.5T solcnoidal magnetic field. The next 
device radially outward is time-of-flight system to distinguish pions, kaons and protons. 
This system only works for lower momenta. The most important advance in the CLEO 
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III, BELLE and BABAR detectors is much better charged hadron identification. Each 
experiment uses different techniques based on Cherenkov radiation to extend ix/K sepa¬ 
ration up to the limit from B decays. The next device is an electromagnetic calorimeter 
that uses Thallium doped Csl crystals; indeed this was the most important new technical 
implementation done in CLEO II and has also been adopted by BABAR and BELLE. 
Afterwards there is segmented iron that serves as both a magnetic flux return and a filter 
for muon identification. 

Figure [21] shows a view of the BELLE detector parallel to the beam. 



Figure 21. through the BELLE detector. 


3.4 b Production Characteristics at Hadron Colliders 

To make precision measurements, large samples of b’s are necessary. Fortunately, these are 
available. With the Fermilab Main Injector, the Tevatron collider will produce « 4 x 10 11 
b hadrons/10 7 s at a luminosity of 2 x 10 32 cm -2 s -1 . These rates compare very favorably to 
e + e~ machines operating at the T(4S). At a luminosity of 10 34 cm _2 s _1 they would produce 
2 x 10 8 B’s/10' s. Furthermore B s , A& and other 6-flavored hadrons are accessible for study 
at hadron colliders. The LHC has about a five times larger b production cross-section. 
Also important are the large charm rates, ~10 times larger than the b rate. 
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In order to understand the detector design it is useful to examine the characteristics 
of b quark production at pp collider. It is often customary to characterize heavy quark 
production in hadron collisions with the two variables pt and 77 , where 77 = —In (tan (6/2 )), 
and 6 is the angle of the particle with respect to the beam direction. According to QCD 
based calculations of b quark production, the B 7 s are produced “uniformly” in 77 and 
have a truncated transverse momentum, p t , spectrum characterized by a mean value 
approximately equal to the B mass (Artuso 1994). The distribution in 77 is shown in 
Figure |22](a). Note that at larger values of 1 77 1, the B boost, /Ty, increases rapidly (b). 

The “flat” p distribution hides an important correlation of bb production at hadronic 
colliders. In Figure ^2|(c) the production angles of the hadron containing the b quark is 
plotted versus the production angle of the hadron containing the b quark according to 
the Pythia generator. Many important measurements require the reconstruction of a b 
decay and the determination of the flavor of the other b, thus requiring both b 7 s to be 
observed in the detector. There is a very strong correlation in the forward (and backward) 
direction: when the B is forward the B is also forward. This correlation is not present in 
the central region (near 90°). By instrumenting a relative small region of angular phase 
space, a large number of bb pairs can be detected. Furthermore the B 7 s populating the 
forward and backward regions have large values of /Ly. 



v 




Figure 22. (a) The B yield versus 77 . (b) /Ly of the B versus p. (c) The production 

angle (in degrees) for the hadron containing a b quark plotted versus the production angle 
for a hadron containing b quark. 


BTeV, a dedicated heavy flavor experiment approved to run at the Fermilab Tevatron 
collider, uses two forward spectrometers (along both the p and p directions) that utilize 
the boost of the B 7 s at large rapidities. This is of crucial importance because the main way 
to distinguish b decays is by the separation of decay vertices from the main interaction. 
LHCb, approved for operation at the LHC, needs a larger detector to analyze the higher 
momentum decay products, and thus has only one arm. 
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3.4.1 The BTeV Detector Description 

I will describe BTeV here though LHCb shares many of the same considerations. There 
are difficulties that heavy quark experiments at hadron colliders must overcome. First of 
all, the huge b rate is accompanied by an even larger rate of uninteresting interactions. At 
the Tevatron the 6-fraction is only 1/500. In searching for rare processes, at the level of 
parts per million, the background from b events is dominant. (Of course all b experiments 
have this problem.) The large data rate of 6’s must be handled. For example, BTeV, 
has 1 kHz of 6’s into the detector, and these events must be selected and written out. 
The electromagnetic calorimeter must be robust enough to deal with the particles from 
the underlying event and still have useful efficiency. Furthermore, radiation damage can 
destroy detector elements. 
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Figure 23. Schematic of the BTeV detector. 

The BTeV Detector is shown in Figure |23] (Skwarnicki 2001) and the LHCb detector 
in Figure ^4| (Muheim 2001). The central part of the BTeV detector has a silicon pixel 
detector inside a 1.5 T dipole magnet. The LHCb experiment uses silicon strips. The 
BTeV pixel detector provides precision space points for use in both the offline analysis and 
the trigger. The pixel geometry is sketched in Figure |25|(a). Pulse heights are measured on 
each pixel. Prototype detectors were tested in a beam at Fermilab; excellent resolutions 
were obtained, especially when reading out pulse heights (Appel 2001) (see Figure |25|(b) ). 
The final design uses a 3-bit ADC. 

The pixel tracker provides excellent vertex resolution, good enough to trigger on events 
with detached vertices characteristic of b or c decays. BTeV shows a rejection of 100:1 for 
minimum bias events in the first trigger level while accepting about 50% of the usable b 
decays. A good explanation of the trigger algorithm can be found at 
http://www-btev.fnal.gov/public_documents/animations/Animated_Trigger^ | . Further 
trigger levels reduce the background by about a factor of twenty while decreasing the b 
sample by only 10%. The trigger system stores data in a pipeline that is long enough to 
ensure no deadtime. The data acquisition system has sufficient throughput to accommo¬ 
date an output of 1 kHz of 6’s, 1 kHz of c’s and 2 kHz of junk. Tracking is accomplished 
using straw tube wire chambers with silicon strip chambers in the high track density 
region near the beam. 

Charged particle identification is done using a Ring Imaging CHerenkov detector. A 
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1) Vertex detector 2) Aerogel and Gas RICH's 3) Magnet yoke 4) Coils 
5) Magnetic field shielding plates 6) Tracking chambers 7) Gas RICH 
8) Electromagnetic calorimeter 9) Hadron calorimeter 10) Muon system 


Figure 24. A schematic diagram of the LHCb detector. 



Figure 25. (a) Pixel detector geometry in BTeV. The detector is inside the beam pipe, 

(b) The spatial resolution as a function of the incident track angle for both 2-bit and 8-bit 
ADC’s as measured in an 800 GeV/c pion beam. The straight lines are piecewise fits to 
the data used in the Monte Carlo simulation. The dashed line near the top indicates the 
resolution obtainable without using pulse height information. 


gaseous C 4 F 10 radiator is used with a large mirror that focuses light on plane of photon 
detectors; these currently are ffybrid Photo-Diodes. They have a photocathode and a 20 
KV potential difference between the photocathode and a silicon diode that is segmented 
into 163 individual pads. The photoelectron is accelerated and focused onto the diode 
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yielding position information for the initial photon. The system will provide four standard- 
deviation kaon/pion separation between 3-70 GeV/c, electron/pion separation up to 22 
GeV/c and pion/muon separation up to 15 GeV/c. Because protons don’t radiate until 9 
GeV/c they can’t be distinguished from kaons below this momenta. BTeV is considering 
an additional liquid CeF 14 radiator, 1 cm thick, in the front of the gas along with a 
proximity focused phototube array adjacent to the sides of the gas volume, to resolve this 
ambiguity. 

In addition, BTeV has an excellent Electromagnetic calorimeter made from Pb\V0 4 
crystals, based on the design of CMS. Finally, the Muon system is used to both identify 
muons and provide an independent trigger on dimuons (BTeV 2000). 


4 b° — B° Mixing 

4.1 Introduction 

Neutral B mesons can transform to their anti-particles before they decay. The diagrams 
for this process are shown in Figure ^6] for the B There is a similar diagram for the B s . 
Although u, c and t quark exchanges are all shown, the t quark plays a dominant role 
mainly due to its mass, as the amplitude of this process is proportional to the mass of 
the exchanged fermion. 



Figure 26. The two diagrams for Bd mixing. 


Under the weak interactions the eigenstates of flavor, degenerate in pure QCD can 
mix. Let the quantum mechanical basis vectors be {11), |2)} = {\B°), |B )}. Then the 
Hamiltonian is 


H = M - l -V 


M Mi 2 
M*2 M 



Diagonalizing we have 


A m = m B[1 - m BL = 2 M u 


(31) 

(32) 


Here H refers to the heavier and L the lighter of the two weak eigenstates. 

Bd mixing was first discovered by the ARGUS experiment (Albrecht 1983) (There was 
a previous measurement by UAl indicating mixing for a mixture of B% and B° (Albajar 
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1987) At the time it was quite a surprise, since m t was thought to be in the 30 GeV 
range. It is usual to define R as probability for a B° to materialize as a B° divided by the 
probability it decays as a B°. The OPAL data for R (Akers 1995) are shown in Figure p7j . 



Figure 27. The ratio, R, of like-sign to total events as a function of proper decay time, 
for selected B D* + X£~P events. The jet charge in the opposite hemisphere is used to 
determine the sign correlation. The curve is the result of a fit to the mixing parameter. 

Data from many experiments has been combined by “The LEP Working Group,” to 
obtain an average value A= (0.489 ± 0.0008) x 10 12 7is -1 . Values from individual 
experiments are listed in Figure 

The probability of mixing is given by (Gaillard 1974) (Bigi 2000) 

^ E y = ^ B Bfl m BT B \V; b V td \ 2 m 2 t F^^^riQCD, (33) 

where Bb is a parameter related to the probability of the d and b quarks forming a hadron 
and must be estimated theoretically, F is a known function which increases approximately 
as m 2 , and Pqcd is a QCD correction, with value about 0.8. By far the largest uncertainty 
arises from the unknown decay constant, /#. In principle fs can be measured. The decay 
rate of the annihilation process B~ —> i~v is proportional to the product of /||I4ft| 2 - Even 
if V u b were well known this is a very difficult process to measure. Our current best hope 
is to rely on unquenched lattice QCD which can use the measurements of the analogous 
D + —> p + v decay as check. This will take the construction of a “r-charm factory.” 

Now we relate the mixing measurement to the CKM parameters. Since 

I V t lV td \ 2 oc |(1 -p- ir/)\ 2 = (p- l) 2 + f] 2 , (34) 

measuring mixing gives a circle centered at (1,0) in the p — rj plane. This could in principle 
be a very powerful constraint. Unfortunately, the parameter Bb is not experimentally 
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ALEPH * 
(3 + 1 prel) 

DELPHI * 
(4 + 1 prel) 


L3 

(3) 

OPAL 

(5) 

SLD * 
(5 prel) 

CDF* 
(4 + 2 prel) 

BABAR * 
(2 prel) 

BELLE 

(1) 


average of above 
after adjustments 

ARGUS+CLEO 
(X d measurements) 

world average 


0.446±0.020±0.018 ps' 1 
0.504±0.019±0.012 ps' 1 
0.444±0.028±0.028 ps' 1 
0.479±0.018±0.015 ps' 1 
0.507±0.023±0.019 ps' 1 
0.495±0.026±0.025 ps' 1 
0.503±0.009±0.011 ps' 1 
0.463±0.008±0.016 ps' 1 

0.48910.009 ps' 1 
0.48910.032 ps' 1 

0.48910.008 ps' 1 

0.4 0.45 0.5 0.55 



working group average .j 

without adjustments Alll^ ) 


Figure 28. Values of the B d mixing parameter A rrtd for each experiment. 


accessible and f B , although in principle measurable, has not been and may not be for 
a very long time, so it too must be calculated. The errors on the calculations are quite 
large. 


4.2 B s Mixing in the Standard Model 


B° mesons can mix in a similar fashion to B% mesons. The diagrams in Figure [26| are 
modified by substituting s quarks for d quarks, thereby changing the relevant CKM matrix 
element from V tc i to V ts . The time dependent mixing fraction is 


Am. s 



gI r 

6n 2 


B, 


f 2 B s m B T Ba \V: b V t 


ts 




Vqcd, 


(35) 


which differs from equation (|33|) by having parameters relevant for the B s rather than the 
B d . 

Measuring x s allows us to use ratio of Xd/x s to provide constraints on the CKM 
parameters p and g. We still obtain a circle in the (p,rj) plane centered at (1,0): 


W, 


td\ 


= A z \ 


2 \ 4 


(i -pY + v 2 


(36) 
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Now however we must calculate only the SU(3) broken ratios B Bd /B Bs and f Bd /f Ba . 

B° mixing has been searched for at LEP, the Tevatron, and the SLC. A combined 
analysis has been performed. The probability, V(t) for a B s to oscillate into a B s is given 
as 



(37) 


where t is the proper time. 

To combine different experiments a framework has been established where each exper¬ 
iment finds a amplitude A for each test frequency u, defined as 



(38) 


Figure ^ shows the world average measured amplitude A as a function of the test fre¬ 
quency u = A m s (Leroy 2001). For each frequency the expected result is either zero for 
no mixing or one for mixing. No other value is physical, although allowing for measure¬ 
ment errors other values are possible. The data do indeed cross one at a A m s of 16 ps -1 , 
however here the error on A is about 0.6, precluding a statistically significant discovery. 
The quoted upper limit at 95% confidence level is 14.6 ps -1 . This is the point where the 
value of A plus 1.645 times the error on A reach one. Also indicated on the figure is the 
point where the error bar is small enough that a 4cr discovery would be possible. This is 
at Am s = 11 ps^ 1 . Also, one should be aware that all the points are strongly correlated. 

The upper limit on A m s translates to an upper limit on x s <21.6 also at 95% confidence 
level. CDF plans to probe higher sensitivity and eventually LHCb and BTeV can reach 
values of 80. 

5 Rare b Decays 

5.1 Introduction 

These processes proceed through higher order weak interactions involving loops, which are 
often called “Penguin” processes, for unscientific reasons (Lingel 1998). A Feynman loop 
diagram is shown in Figure ^U| that describes the transition of a b quark into a charged 
-1/3 s or d quark, which is effectively a neutral current transition. The dominant charged 
current decays change the b quark into a charged +2/3 quark, either c or u. 

The intermediate quark inside the loop can be any charge +2/3 quark. The relative 
size of the different contributions arises from different quark masses and CKM elements. 
In terms of the Cabibbo angle (A=0.22), we have for t:c:u - A 2 : A 2 :A 4 . The mass dependence 
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Figure 29. Combined experimental values of the amplitude A versus the test frequency 
to = A rrid as defined in equation |i?i|. The inner (outer) envelopes give the 95% confidence 
levels using statistical (statistical and systematic) errors. The “sensitivity” shown at 18.3 
ps^ 1 is the likely place a 95% c.l. upper limit could be set. Also indicated is the maximum 
value, 11 ps~ l , where a fo discovery would be possible. 


W" 



Figure 30. Loop or “Penguin” diagram for a b —> s orb^d transition. 

favors the t loop, but the amplitude for c processes can be quite large ~30%. Moreover, as 
pointed out by Bander, Silverman and Soni (Bander 1979), interference can occur between 
t, c and u diagrams and lead to CP violation. In the Standard Model it is not expected 
to occur when b —> s, due to the lack of a CKM phase difference, but could occur when 
b —> d. In any case, it is always worth looking for this effect; all that needs to be done, 
for example, is to compare the number of K*~ 7 events with the number of K* + 7 events. 

There are other ways for physics beyond the Standard Model to appear. For example, 
the W~ in the loop can be replaced by some other charged object such as a Higgs; it is 
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also possible for a new object to replace the t. 

5.2 Standard Model Theory 

In the Standard Model the effective Hamiltonian for the intermediate t quark is given by 
(Desphande 1994) 

4C 10 

He,f = - JVttl. (39) 

V 4 i =i 

Some of the more important operators are 

Oi = s\rf,MAYc J L. 0 7 = m h s l L a, w I jP jF Nl . (40) 

The matrix elements are evaluated at the scale /i = M w and then evolved to the b 
mass scale using renormalization group equations, which mixes the operators: 


C i (fi) = J2U l MM w )C j (M w ) . (41) 

3 


5.3 b —► sy 

This process occurs when any of the charged particles in Figure [30] emits a photon. The 
only operator which enters into the calculation is C 7 (fi). We have for the inclusive decay 


Heff 


O 


7 


T (6 —► s'y) 


^=-(V tb V ts \ 2 C 7 (m b )0 7 

(42) 

16?r 

(43) 

G 2 F am 5 b 2 2 

327T4 |C71 lVtbVt ° l ■ 

(44) 


Its far more difficult to calculate the exclusive radiative decay rates, but they are much 
easier to measure. Note that the reaction B —> K'y would violate angular momentum 
conservation, so the simplest exclusive final states are B —» K* 7 . 

Different techniques are used for reconstructing exclusive and inclusive Decays and 
unique methods are invoked for exclusive decays on the T(4S'). At other machines the 
decay products, i, from an exclusive B decays are used to reconstruct an “invariant 
mass” via M 2 = At the Y(4S) its done a bit differently, the decay 

products are first tested to see if the sum of their energies is close to the beam energy, 
Ebeam■ If this is true then the “beam constrained” invariant mass is calculated as M 2 = 
Eb eam ~ Si 'P’i ■ These methods are used for all exclusive B decays, in combination with 
other requirements. Figure [3T] shows the BELLE data for the reaction B D* + n~, 
where the D* + —> tt + D°. BELLE Erst selects events with candidate D°' s. Then they 
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require an additional 7 r + where the measured mass difference between the D°tt + minus the 
D° candidate is consistent with the known mass difference. Selecting the D* + candidates 
they combine them with candidate n~. In Figure [H](a) the correlation between difference 
in measured energy A E = Yii E t — Eb eam versus the beam-constrained invariant mass 
is shown. In Figure [JT](b) A E is shown after selecting the signal region in M and in 
Figure |3l|(c) M is shown after selecting on A E. These plots show how clean signals can 
be selected. 


B°-> D V 



Figure 31. BELLE data for the reaction B° —> D* + 7r~. (a) The correlation between A E 
and M. The box shows the signal region, (b) The projection in A E for events in the M 
signal region; the line shows a fit to the background, (c) The M distribution for A E in 
the signal region; the line shows a fit to the background. 

CLEO first measured the inclusive rate (Alam 1995) as well as the exclusive rate into 
K* (890)7 (Amniar 1993) shown in Figure [32]. Here several different decay modes of the 
K*{ 890) are used. The current world average value for B(B —* K*y) = (4.2 ±0.8 ±0.6) x 
10 ~ 5 . 

To End inclusive decays two techniques are used. The first one, which provides the 
cleanest signals, is to sum the exclusive decays for the final states Krnry, where n < 4 and 
only one of the pions is a 7T°. These requirements are necessary or the backgrounds become 
extremely large. (Both charged and neutral kaons are used.) Of course, imposing these 
restrictions leads to a model dependence of the result that must be carefully evaluated. 
This is why having an independent technique is useful. That is provided by detecting only 
the high energy photon. The technique used is to form a neutral network to discriminate 
between continuum and T(4S') data using shape variables. 

The momentum spectrum of the 7 peaks close to its maximum value at half the B 
mass. If we had data with only B mesons, it would be easy to pick out b —► sy. We have, 
however, a large background from other processes. At the T(4S), the 7 spectrum from 
the different background processes is shown. The largest is n° production from continuum 
e + e _ collisions, but another large source is initial state radiation (ISR), where one of the 
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5,200 5.220 5,240 5.260 5,280 5.300 

U (GeV) 

Figure 32. First published CLEO data for the reaction B —> K *7 showing the M 
distribution for A E in the signal region. 


beam electrons radiates a hard photon before annihilation. The backgrounds and the 
expected signal are illustrated in Figure |33|. Similar backgrounds exist at LEP. 



Figure 33. Levels of inclusive photons from various background processes at the T(4S') 
labeled largest to smallest at 2.5 GeV/c. Also shown is the expected signal from b —> s'y. 

To remove background CLEO used two techniques originally, one based on “event 
shapes” and the other on summing exclusively reconstructed B samples. Examples of 
idealized events are shown in Figure [34|. CLEO uses eight different shape variables de¬ 
scribed in Ref. [3], and defines a variable r using a neural network to distinguish signal 
from background. The idea of the B reconstruction analysis is to find the inclusive branch- 
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Continuum Continuum +ISR Signal 



Figure 34. Examples of idealized event shapes. The straight lines indicate hadrons and 
the wavy lines photons. 


ing ratio by summing over exclusive modes. The allowed hadronic system is comprised 
of either a K s —> 7r + 7r~ candidate or a K T combined with 1-4 pions, only one of which 
can be neutral. The restriction on the number and kind of pions maximizes efficiency 
while minimizing background. It does however lead to a model dependent error. For all 
combinations CLEO evaluates 


Xb = 


M, 


B 


&M 


5.279 y 


+ 


E r — E, 


beam 


(JE 


; ■ 


(45) 


where M B is the measured B mass for that hypothesis and E B is its energy. Xb i s 
required to be < 20. If any particular event has more than one hypothesis, the solution 
which minimizes Xb is chosen. For events with a reconstructed B candidate CLEO also 
considers the angle between the direction of the B and the thrust axis of event with the B 
candidate removed, cos(6^). This is highly effective in removing continuum background. 


A neural network is used to combine r, Xbi cos (9 t ) into a new variable r c and events 
are then weighted according to their value of r c . This method maximizes the statistical 
potential of the data. Figure 35 shows the photon energy spectrum of the inclusive 
signal from CLEO combining both reconstruction techniques. The signal is compared to 
a theoretical prediction based on the model of Ali and Greub (Ali 1991). A fit to the 
model over the photon energy range from 2.0 to 2.7 GeV/c gives the branching ratio result 
shown in Table [5|, where the first error is statistical, the second systematic and the third 
dependence on the theoretical model (Chen 2001). 


ALEPH reduces the backgrounds by weighting candidate decay tracks in a 6 —> sy 
event by a combination of their momentum, impact parameter with respect to the main 
vertex and rapidity with respect to the b- hadron direction (Barate 1998). 

Current results are also shown in Table [5|. The data are in agreement with the Standard 
Model theoretical prediction to next to leading order, including quark mass effects of 
(3.73 ± 0.30) x 10~ 4 (Hurth 2001). A deviation here would show physics beyond the 
Standard Model. More precise data and better theory are needed to further limit the 
parameter space of new physics models, or show an effect. 
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1.5 2.5 3.5 4.5 


E y (GeV) 

Figure 35. The background subtracted photon energy spectrum from CLEO. The spectrum 
is not corrected for resolution or efficiency. The solid lines show the spectrum from a 
simulation of the Ali-Greub spectator model with the b quark mass set to 4.690 GeV and 
the Fermi momentum set to 410 MeV/c. 


Table 5. B{b —> s'y). 


Experiment 

B x 10~ 4 

CLEO 

3.21 ±0.43 ±0.27ig;lg 

ALEPH 

3.11 ±0.80 ±0.72 

BELLE 

3.36 ± 0.53 ± 0.44^ 

Average 

3.23 ±0.42 


5.3.1 \Vcb\ Using Moments of the Photon Energy Spectrum 


In section (|2.2.4| ) we found a value of V c b using the first and second hadronic mass moments. 
Here we use the first moment of the photon energy distribution in b —> s'y. The values 
found for the moments and A which is directly proportional to (Ef) are (Chen 2001) 


(EJ 

(E*) ~ (E,f 

A 


2.346 ±0.032 ±0.011 GeV 
0.0226 ± 0.0066 ± 0.0020 GeV 2 
0.35 ±0.08 ±0.10 GeV. 


(46) 

(47) 

(48) 


In Figure ^ we show the combination of first moments from photon energy in b —> s'y 
and hadron moments in b —> civ. This implies a value of V c b around 0.0406. 
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- 0.4 

- 0.5 

A 

Figure 36. Correlation between Ai, A and V c b derived from (Ef) and — M D ). The 
lighter bounds include both experimental and theoretical errors. 



5.4 Rare Hadronic Decays 
5.4.1 B — » 7T7T and B — > K tv 

The decays B° —> 7t + tt~ and B° — > K~tt + do not contain any charm quarks in the final 
state so must proceed via either the tree level V u b process shown in Figure ^7|(left) or via 
the Penguin process shown on the right side. 



(right) Via a Penguin process. 


w 



(left) Via tree level V u b moderated decay. 


These diagrams can interfere and they can also interfere through B° mixing, thus 
complicating any weak phase extraction. The same diagrams are applicable for B —> 
K~7 r + by replacing W~ —> ud in the tree level diagram by W~ —> us and replacing the 
td coupling in the Penguin by a ts coupling. 


Other diagrams for producing Krc final states are shown in Figure 


In section 7.1 
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it will be shown that CP violation can result from the interference between two distinct 
decay amplitudes leading to the same final state. Consider the possibility of observing 
CP violation by measuring a rate difference between B~ —> K~tt° and B + —> K + 7r°. The 
K~7i° final state can be reached either by tree or penguin diagrams. The tree diagram 



Figure 38. Diagrams for B~ —> K~ir°, (a) and (b) are tree level diagrams where (b) 
is color suppressed; (c) is a penguin diagram, (d) shows B~ —> K°tt~ , which cannot be 
produced via a tree diagram. 

has an imaginary part coming from the V u b coupling, while the penguin term does not, 
thus insuring a weak phase difference. This type of CP violation is called “direct.” Note 
also that the process B~ —> K°tt~ can only be produced by the penguin diagram in 
Figure ^(d). Therefore, we do not expect a rate difference between B~ —> K°tt~ and 
B + -+ K°tt + . 

Measurements of these rates have been by several groups. Recent data from BELLE 
are shown in Figure |3^ (Abe 2001a). Table lists the currently measured branching 
ratios. 


Table 6. Branching Ratios for B —> Kn and B —> mr in units of 10 


Mode 

CLEO 

BABAR 

BELLE 

Average 

7T + 7T - 

4.7ll;| ±0.6 

4.1 ±1.0 ±0.7 

5.6ll:g ±0.4 

4-5+1 

7r + 7T° 

<12 

<9.6 

<13.4 


K ±r K T 

18.811:1 ± 1-3 

16.7 ± 1.6 ± 1.3 

i q o+3.4+1.5 
3.2—.06 

17.71+ 

K + tt° 

i o 1 +3.0+2.1 

1^. 1—2.8—1.4 

10.8+9 ± 1-6 

1 f? q+3.5+1.6 
i0 -°-3.3-1.8 

12.H+ 

o 

n 

i 

i8.2lt;B ± 1.6 

18.2ll;l ±2.0 

i O 7+5. 7+1.9 
' -4.8-1.8 

17.311:1 

K°tt° 

i4.8l| : ?li:t 

8. 2!+ ± 1.2 

i0 -U 5.9 2 .7 

10.41+ 
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M (GeV/c 2 ) AE (GeV) M (GeV/c 2 ) AE (GeV) 


Figure 39. Signals in M and AE for two-body decay modes from BELLE. The data 
result from projections of a Likelihood fit that takes into account event shape and particle 
identification information. The dashed lines are the signal projections. The dotted lines 
in the AE distributions are projections of the background component from the n K 
substitution. 

6 Hadronic Decays 

6.1 Introduction 

Mark Wise in his talk at the 2001 Lepton Photon conference (Wise 2001) gave some 
advice to theorists: “If you drink the nonlep tonic your physics career will be ruined and 
you will end up face down and in the gutter.” Presumably Mark’s statement was inspired 
by the difficulty in predicting hadronic decays. Here we have lots of gluon exchange 
with low energy gluons, while perturbation theory works well when the energies are large 
compared with A qcd ~200 MeV. Furthermore, multibody decays are currently impossible 
to predict, so we will consider only two-body decays. 

6.2 Two-Body Decays into a Charm or Charmonium 

We start by considering two-body decays into a charmed hadron (Neubert 1997). Figure 
shows the processes for both B~ and B decays into a D and a . There is only one 
possible process for the B , the simple spectator process (left), while the color suppressed 
spectator (right) is also allowed for the B~. We call decays with only the simple spectator 
diagram allowed “class I” and decays where both the simple and color suppressed diagrams 
are allowed “class III.” Note, that because the colors of all the outgoing quarks must be the 
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same in the color suppressed case, naively the amplitude is only -—'1/3 that of the simple 
spectator case where the W~ can transform into quarks of all three colors. “Class II” 
decays are processes than can only be reached by the color suppressed spectator diagram, 
for example the B c 


J / if K s decay shown in Figure [40 



Figure 40. Color suppressed spectator decay diagram for B° —> J/ifK s . 


The effective Hamiltonian consists of local 4-quark operators renormalized at the scale 
H and the Wilson coefficients (from the Operator Product Expansion) We have 


Heff 

Qi 

Qt 


cb c \^h)Q c \ S-Cip)QiCb + Wb---} 

(d u)v-A + (s'c)v-A ( c'b)v-A 
C cu) V - A (d'b) V -A + ( cc)v-A(s'b)v_A , 


(49) 


where the notation (' q^q 2 )v-A = gT7 M (l — 75 )g 2 - Without QCD corrections c\(fi) = 1 and 
c 2 (/u) = 0- For non-leading order correction using the renormalization group equations, 
we have Ci(/i) = 1.132 and c 2 (/i) = —0.249. 

We can factorize the amplitude by assuming that the current producing the n~ is 
independent of the one producing the charmed hadron. Lets consider a class I case, 
B D + tt~ . The amplitude can be written as 


A i^ = v a v:m* 


(du) A 0)(D+\(cb) v \B°) 


(50) 


The part of the amplitude dealing with the n is known from pion decay. We have 
(tt~ d'y^u 0) = ifnPfj,, where the axial vector structure is made explict, p M is the pion 
four-vector and f n is given by measuring the decay width for n~ —► p~v. The term a\ is 
defined as 

ai = c 1 {nf)+£c 2 {nf) , (51) 


where £ is equal to the number of colors and the scale pj is on the order of the b quark 
mass. Then 

A l~* = I - ml)F 0 B ^ D (ml) . (52) 

The Fq form factor can either be calculated or measured, for example in semileptonic 
decays. 
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Let us also consider a class II process B —> J/ipK. In this case 

A fact = -^Lv cb V*a 2 {J/i;\(cc)v\ 0)(K\(sb) v \B) , (53) 

where a 2 = c 2 (/i/) + £ci(/i/). 

A class III decay has a term oq + xa 2 in the amplitude, where x equals one from flavor 
symmetry. However the actual values of ai, x, and a 2 are not well predicted from theory 
but we can obtain them from the data. 

One method is to use the class I decays to obtain cq = 1.08 ± 0.04. It is possible to 
calculate a 2 /a\ as shown in Figure [|l] (Neubert 1997). Using these values the measured 
branching ratios are compared with the predicted ones in Table |7]. Here x = +1 is used, 
taken from the data. This is opposite to the interference in D decays but is expected from 
the calculation shown in Figure |4l]. 



C^sCM-f) 

Figure 41. Calculated values of a 2 /a\ versus a s . The arrows indicate the ratios chosen 
for B and D decays. 

Table 7. Predicted Branching Ratios (T0 -4 j (Neubert 1997) Compared to Measurement 
For Two-Body Hadronic Decays. 



Class I 



Class II 



Class III 

Mode 

Model 

Data 

Mode 

Model 

Data 

Mode 

Model 

Data 

D + 7T- 

30 

30±4 

J/ipK- 

11 

10±1 

D°tt~ 

48 

53±5 

D* + 7T~ 

30 

28±2 

J/ifK* 

17 

15±2 

D*°7T~ 

49 

46±4 

D+p- 

70 

79±14 

D°tt° 

0.7 

2.8T0.4 

D°p~ 

110 

134±18 

D* + p- 

85 

73±15 

D*°tt° 

1.0 

2.1T0.9 

D*°p~ 

119 

155±31 


The agreement is rather good except for the newly measured D°^*\° modes where it 
is rather miserable. CLEO and BELLE both have rates for D°n° of (2.7 ± 0.3 ± 0.5) x 
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10 4 and (2.91°;! ± 0.6) x 10 4 , respectively, while CLEO alone has measured D*°ir° as 
(2.1 ±0.5 ±0.8) x 10~ 4 . 

6.2.1 Isospin Analysis of the B —> Dn System 

All the decay rates for B —> Dix have now been measured. The four-quark operator 
( du)(cb ) has isospin 1=1 and l 3 =+l. It transforms the B into final states D + n~ and 
D°ir° with 1=1/2 or 1=3/2. The B~ decays into D°tt~ with 1=3/2 only. It is thought 
that the isospin amplitudes cannot be modified by final state interactions, so we can 
look for evidence of final state phase shifts by doing an “isospin analysis.” The decay 
amplitudes form a triangle as shown in Figure f|^. 



Figure 42. The isospin amplitude triangle for B —> Dtt decays. 

The relationship among B decay amplitudes and isospin amplitudes is given by 


A(B° -> D + n~) 

A(B° -»• D°tt°) 
A(B~ -»■ D°7r~) 



(54) 


These equations may be solved for the isopsin amplitudes and the relative phase shift 
between the two amplitudes. The solution is 


A 


1/2 


A 


1/2 


A(B° D + tt~ 


A(B~ -> D°ti~) 


A(B° 


D°tt°)\ 2 - ^\A(B- 


D°ir~ 


(55) 


, \ 3\A(B~ -> D°7t-)| 2 -2 

A. ,A — 1 

-^1/2 

2 

^3/2 

CO 

tc 

C 

1 

bO 

h, 

x4-1/2 

^3/2 



Solving these equations for the Dtx final states gives cosh = 0.88±0.05 which indicates 
a phase shift of about 28 ± 8 degrees, but is not statistically significant. 
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6.2.2 Factorization Tests Using Semileptonic Decays 


The factorized amplitude for Dn~ decays in equation [53] is the product of two hadronic 
currents, one for W~ and the other for B —> D. In semileptonic decay (Figure ||]) 

we have the product of the known lepton current and the pion current. At q 2 = m 2 the 
B —> D should be the same in both decays, at least for class I. The comparision for the 
general case of any hadron h~ is 


r (B -»■ D^h) 


6 K 2 a\fl \V ud \ 




D^n 7 ) 


q 2 =m% 


(56) 


Tests of this equation for D* + and a 7 r , p or a 1 are satisfied at about 15% accuracy 
(Bortoletto 1990) (Browder 1996). 


Another test compares the polarization of the D* in both hadronic and semileptonic 
cases: 


Tl 

r 









(57) 


where T L denotes the longitudinally polarized fraction of the decay width. Comparisons 
with data will be shown in section |6.3| . 

There are also more modern approaches to factorization (Beneke 2001) (Bauer 2001). 
However these approaches predict 


r (b~ -> d°tt-) 
r (b° -> d+tt-) 


1 + 0(\QCD/n^b ) 


which seems to contradict current observations. 


(58) 


6.3 Observation of the p f in B Decays 

CLEO made the first statistically significant observations of six hadronic B decays shown 
in Table |8| that result from studying the reactions B —> (Alexander 2001). 

The signal in one of these final states B —t ZA* + )7r + 7r _ 7r _ 7r°, where D* + —> tt + D° and 
D° —> K~7i + is shown in Figure [43], 

In examining the substructure of the four-pions, a clear u signal was observed in the 
n + 7i~ir° mass as can be seen in Figure [0], leading to a significant amount of t~. 

Furthermore, there is a low-mass resonant substructure in the ujh~ mass. (See Figure f45[). 

The spin and parity of the iV7i~ resonance (denoted by A temporarily) is determined 
by considering the decay sequence B —> A D] A —> un and ui —> 7r + 7r _ 7r°. The angular 
distributions are shown in Figure £h| Here 6a is the angle between the u direction in the 
A rest frame and the A direction in the B rest frame; 9 b is the orientation of the u decay 
plane in the u rest frame, and y is the angle between the A and u decay planes. 
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Figure 43. The B candidate mass spectra for the final state D* +/ k + ti~ti~'k 0 , with D° —> 
K~ 7r + (a) for A E sidebands arid (b) for A E consistent with zero. The curve in (a) is a 
fit to the background distribution described in the text, while in (b) the shape from (a) is 
used with the normalization allowed to float and a signal Gaussian of width 2.1 MeV is 
added. 


Table 8. Measured Branching Ratios 


Mode 


B (%) 

B° 

D* + 7r + 7r _ 7r _ 7r° 

1.72±0.14±0.24 

B° -»■ 

D* + C07l~ 

0.29±0.03±0.04 

B° -> 

D + ujtt~ 

0.28±0.05±0.03 

B~ -> 

■ D*°Tr + 7T~7T~7r 0 

1.80±0.24±0.25 

B~ -» 

■ D*°um~ 

0.45±0.10±0.07 

B~ -» 

■ D°um~ 

0.41±0.07±0.04 


The data are fit to the expectations for the various J p assignments. The u> polarization 
is very clearly transverse (sin 2 9 U ) and that infers a 1“ or 2 + assignment. The cos 9a 
distribution prefers 1 _ , as does the fit to all three projections. 

This structure is identified with the p' because it has the correct J p and is at approx- 
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Figure 44. The invariant mass spectra o/7r + 7r _ 7r° for the final state D* + tt + 7i~'ti~7i 0 for 
all three D° decay modes K~n + , K~7r + n + n~ and K~tt + tt 0 . The solid histogram is the 
background estimate from the Mb lower sideband and the dashed histogram is from the 
A E sidebands; both are normalized to the fitted number of background events. (There is 
an additional cut selecting the center of the Dalitz plot.) 


1630301-037 



Figure 45. (left) The background subtracted efficiency corrected lvtt mass spectrum from 
B° —> ( D* + + D° + Z} + )u;7r - decays fit to a Breit-Wigner shape. 


imately the right mass. To determine the mass and width parameters, that are not well 
known, we write the decay width as a fu n ction of tn7r~ mass as 


dT (B —> Dlott ) = 


'A(B —> Dp') ■ 13W(m U7r ) ■ A(p' —> inn )' 


2 M; 


B 


(59) 


xdP(B -»■ Dp) ■ dP(p -> unr~) 


dm 


2 

LOTT 


2vr 
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Figure 46. The angular distribution of 6 a (top-left),9 U (top-right) and x (bottom). The 
curves show the best fits to the data for for different J p assignments. The 0“ and 1 + are 
almost indistinguishable in cos 9a, while the 1“ and 2 + are indistinguishable in cos 9^ and 
X■ The vertical axis gives efficiency corrected events, 10f events are used. 


where dP indicates phase space and the Breit-Wigner amplitude is given by 


6W(m„) 


1 

(mfjix - m 2 pl ) - 


(60) 


The Breit-Wigner fit assuming a single resonance and no background gives a mass of 
1349±25t5° MeV with an intrinsic width of 547±86±45 MeV. The fit shows that the u;7r _ 
mass spectrum is consistent with being entirely one resonance. This state is likely to be 
the elusive p' resonance (Clegg 1994). These are by far the most accurate and least model 
dependent measurements of the p' parameters. The p' dominates the final state. (Thus 
the branching ratios for the apply also for D^p'~.) 

Heavy quark symmetry predicts equal partial widths for D*p' and Dp'. CLEO measures 
the relative rates to be T (B —> D*p'~)/T (B —»• Dp'~) = 1.06 ± 0.17 ± 0.04, consistent 
within the relatively large errors. 

Factorization predicts that the fraction of longitudinal polarization of the D* + is the 
same as in the related semileptonic decay B —> D*l~u at four-momentum transfer q 2 
equal to the mass-squared of the p' 


r L (b D* + p'~) 
F(B -> D*+p'-) 


r L (B -> D*f p) 
T (B -»• D*t-v) 


=m 


2 

p' 


(61) 


CLEO’s measurement of the D* + polarization is (63±9)%. The model predictions in 
semileptonic decays for a q 2 of 2 GeV 2 , are between 66.9 and 72.6% (Isgur 1990) (Wirbel 
1985) (Neubert 1991). Figure [47| shows the measured polarizations for the D* + p '~, the 
D* + p ~, (Artuso 1999) and the D* + D*~ final states (Stone 2000). The latter is based on 
a new measurement using partial reconstruction of the D* + (Ahmed 2000). Thus this 
prediction of factorization is satisfied. 
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Figure 47. Measured D* + polarization versus semileptonic model predictions. 


7 CP Violation 

7.1 Introduction 

Recall that the operation of Charge Conjugation (C) takes particle to anti-particle and 
Parity (P) takes a vector to —~r f . CP violation can occur because of the imaginary 
term in the CKM matrix, proportional to rj in the Wolfcnstein representation (Bigi 2000). 

Consider the case of a process B —* / that goes via two amplitudes A and B each of 
which has a strong part e. g. s .4 and a weak part w_a. Then we have 


r (B 


T(B - 

-/) 

= ( ^4 e i(sA+WA) + \B\ e ^ SB+WB) ) 2 

(62) 

t(b - 

-7) 

= e i{sA ~ WA) + \B\ e ^ SB ~ WB) Y 

(63) 

f ) - r (b - 

-7) 

= 2 \ AB\ sin(s^ — sb) sin(u’_4 — wb ) • 

(64) 


Any two amplitudes will do, though its better that they be of approximately equal 
size. Thus charged B decays can exhibit CP violation as well as neutral B decays. In 
some cases, we will see that it is possible to guarantee that |sin(s _4 — sg)| is unity, so we 
can get information on the weak phases. In the case of neutral B decays, mixing can be 
the second amplitude. 
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7.2 Unitarity Triangles 


The unitarity of the CKM matrixQ allows us to construct six relationships. The most 
useful turns out to be 

VudVtd + V w V{ a + KtUft = 0 . (65) 

To a good approximation 


then 


Since V us 


V ud ~ V* b « 1 and V* s « -Kb, 


Ydh - y = o 
Kb + v* 

A, we can dehne a triangle with sides 


1 

V td 

AX 3 

Vub_ 

AX 3 


f+x 


i 

A 


Kb 

Kb 


Vtd 

V ts 


( 66 ) 

(67) 

( 68 ) 

(69) 

(70) 


The CKM triangle is depicted in Figure |48|. We know the length of two sides already: 



Figure 48. The unitarity triangle shown in the p — r] plane. The left side is determined 
by measurements ofb^ujb^c and the right side can be determined using mixing mea¬ 
surements in the B s and B^ systems. The angles can be found by making measurements 
of CP violating asymmetries in hadronic B decays. 


the base is defined as unity and the left side is determined by the measurements of 
|Kb/Kb|, but the error is still quite substantial. The right side can be determined using 
mixing measurements in the neutral B systems. Figure |4^ also shows the angles as a, f3, 
and 7 . These angles can be determined by measuring CP violation in the B system. 


Wnitarity implies that any pair of rows or columns are orthogonal. 
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Another constraint on p and r; is given by the Kf CP violation measurement (e) (Buras 
1995): 

77 [(1 - p)A 2 {IA ± 0.2) + 0.351 A 2 — A- = (0.30 ± 0.06), (71) 


where Bk is parameter that cannot be measured and thus must be calculated. A reason¬ 
able range is 0.9 > > 0.6, given by an assortment of theoretical calculations (Buras 

1995); this number is one of the largest sources of uncertainty. Other constraints come 
from current measurements on V u b/V c b , Bd mixing and B s mixing. The current status of 
constraints on p and r\ is shown in Figure |49] (Hocker 2001). The width of both of these 
bands are also dominated by theoretical errors. This shows that the data are consistent 
with the standard model. 



Figure 49. The regions inp—rj space (shaded), where p = p( 1—A 2 /2) andrj = rj( 1—A 2 /2), 
consistent with measurements of CP violation in Kf decay ( e ), V u b/V c b in semileptonic B 
decay, B% mixing, and the excluded region from limits on B° mixing. The allowed region 
is defined by a fit from (Hocker 2001). The large width of the Bd mixing band is dominated 
by the uncertainty in .Bb/1- The lines that are not specified are at 5% confidence level. 

ft is crucial to check if measurements of the sides and angles are consistent, i.e., 
whether or not they actually form a triangle. The standard model is incomplete. It has 
many parameters including the four CKM numbers, six quark masses, gauge boson masses 
and coupling constants. Perhaps measurements of the angles and sides of the unitarity 
triangle will bring us beyond the standard model and show us how these paramenters are 
related, or what is missing. 

Furthermore, new physics can also be observed by measuring branching ratios which 
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violate standard model predictions. Especially important are “rare decay,” processes such 
as B —> or D —> 7r/i + /x _ . These processes occur only through loops, and are an 

important class of Penguin decays. 


7.3 Formalism of CP Violation in Neutral B Decays 

Consider the operations of Charge Conjugation, C, and Parity, P: 

C\B(t)) = C\B{?)) = (72) 

P\B{?)) = -| B(-t)), P\B{?)) = -| B{-?)) (73) 

CP\B(?)) = -\B(-?)), CP\B(t)) = -\B{-?)) ■ (74) 

For neutral mesons we can construct the CP eigenstates 

I B{) = 2= (|B°) - |B°>) , (75) 

\B1) = 4= (|B°) + |B°)) , (76) 

where 

CP\B°) = | B{) , (77) 

CP\B°) = -| B°) . (78) 

Since B° and B° can mix, the mass eigenstates are a superposition of a\B°) +b\B°) which 
obey the Schrodinger equation 


.d 

i— 

dt 






(79) 


If CP is not conserved then the eigenvectors, the mass eigenstates | Bj) and \B H ) 1 are not 
the CP eigenstates but are 


Bl) = P \B°) + q\B °), | B h ) = p\B°) - q\B °), 


where 


1 1 + ee ^ 1 1 — 6b 

TSj/rTM 5 ’ q ~ Tlj/TT^P 


CP is violated if e B ^ 0, which occurs if \q/p\ ^ 1. 
The time dependence of the mass eigenstates is 


| B L (t)) = e _rt,/2 e i ’" 1,/2 |B 1 (0)) 
B„(t)) = e- r *‘ /2 e i ’"»*'' 2 |B H (0)), 


(80) 

(81) 


(82) 

(83) 
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leading to the time evolution of the flavor eigenstates as 



(84) 

(85) 


where m = ( m,L + m#)/2, Am = m# — and Y = Yl ~ Yh, and t is the decay time in 
the B° rest frame, the so called proper time. Note that the probability of a B° decay as a 
function of t is given by {B°(t)\B°(t))*, and is a pure exponential, e~ 14 / 2 , in the absence 
of CP violation. 

7.3.1 CP violation for B via interference of mixing and decays 

Here we choose a final state / which is accessible to both B° and B° decays (Bigi 2000). 
The second amplitude necessary for interference is provided by mixing. Figure |50] shows 
the decay into / either directly or indirectly via mixing. It is necessary only that / be 



/ 


Figure 50. Two interfering ways for a B° to decay into a final state f. 

accessible directly from either state; however if / is a CP eigenstate the situation is far 
simpler. For CP eigenstates 


CP\fcp) = ±1 fcp)- 


( 86 ) 


It is useful to define the amplitudes 


A = (f C p\H\B°), A — {fcp\H\B°). 


(87) 


If ^ 1, then we have “direct” CP violation in the decay amplitude, which we will 
discuss in detail later. Here CP can be violated by having 



( 88 ) 


which requires only that A acquire a non-zero phase, i.e. |A| could be unity and CP 
violation can occur. 
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A comment on neutral B production at e + e _ colliders is in order. At the T(4S') 
resonance there is coherent production of B°B° pairs. This puts the B's in a C = — 1 
state. In hadron colliders, or at e + e _ machines operating at the Z°, the R’s are produced 
incoherently. For the rest of this article I will assume incoherent production except where 
explicitly noted. 

The asymmetry, in this case, is defined as 

r -> fcp ) - r (B°(t) -> fcp ) 
a,CF ~ r (s»(t) -> f C p) + r - f C p) ' 

which for \q/p\ = 1 gives 

(1 — |A| 2 ) cos (Amt) — 21mA sin(Amt) 

afcp = l+W 

For the cases where there is only one decay amplitude A, |A| equals 1, and we have 


(89) 


(90) 


df CP = —ImAsin(A mt). 


(91) 


Only the amplitude, —ImA contains information about the level of CP violation, the sine 
term is determined only by Bd mixing. In fact, the time integrated asymmetry is given 
by 

af CP = - --ImA = — 0.481mA . (92) 

1 I oc 

This is quite lucky as the maximum size of the coefficient for any x is —0.5. 

Let us now find out how ImA relates to the CKM parameters. Recall A = | ■ (f- The 
first term is the part that conies from mixing: 


q (V t lV td f (1 ~ P ~ irif 

P \V tb V td \ 2 (1 -p + vri) (1 -p-irj) 



2(1 ~P)V 
(1 ~ P? + V 2 


sin (2 (3). 


(93) 

(94) 


To evaluate the decay part we need to consider specific final states. For example, 
consider / = 7r + 7r~. The simple spectator decay diagram is shown in Figure 37 (left). 
For the moment we will assume that this is the only diagram which contributes, which 
we know is not true. For this b —> uud process we have 


A = (VSjVW 2 = (p - i >)) 2 = e _ 2h 

A \VudVubf (P-»’!)(P + Mf) 

Im(A) = Im(e~ 2i/3 e~ 2 * 7 ) = Im(e 2i “) = -sin(2a). 


(95) 

(96) 
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For our next example let’s consider the final state J/i[)K s . The decay diagram is shown 
in Figure |I]. In this case we do not get a phase from the decay part because 


A (VdK'f 

A IKWcl 2 * 


(97) 


is real to order 1/A 4 . 

In this case the final state is a state of negative CP, i.e. CP\J/iftK s ) = —| J/ipK s ). 
This introduces an additional minus sign in the result for ImA. Before finishing discussion 
of this final state we need to consider in more detail the presence of the K s in the final 
state. Since neutral kaons can mix, we pick up another mixing phase (similar diagrams 
as for B°, see Figure |26|). This term creates a phase given by 


M (v; d v„f 

\P/ k 'V'4 2 


(98) 


which is real to order A 4 . It necessary to include this term, however, since there are other 
formulations of the CKM matrix than Wolfenstein, which have the phase in a different 
location. It is important that the physics predictions not depend on the CKM convention.fi 

In summary, for the case of / = J/x/iK s , ImA = — sin(2/3). 


7.3.2 Comment on Penguin Amplitude 

In principle all processes can have penguin components. One such diagram is shown in 
Figure |37](right). The 7r + 7r - final state is expected to have a rather large penguin ampli¬ 
tude ~10% of the tree amplitude. Then |A| ^ 1 and a 7r7r (f) develops a cos(A mt) term. It 
turns out that sin(2a) can be extracted using isospin considerations and measurements 
of the branching ratios for B + —> 7T + 7r° and B° —> i\°i r°, or other methods the easiest of 
which appears to be the study of B° —> pir. 

In the J/'ipKg case, the penguin amplitude is expected to be small since a cc pair must 
be “popped” from the vacuum. Even if the penguin decay amplitude were of significant 
size, the decay phase is the same as the tree level process, namely zero. 


7.4 Results on sin 2/3 

For years observation of large CP violation in the B system was considered to be one 
of the corner stone predictions of the Standard Model. Yet it took a very long time to 
come up with definitive evidence. The first statistically significant measurements of CP 
violation in the B system were made recently by BABAR and BELLE. This enormous 

2 Here we don’t include CP violation in the neutral kaon since it is much smaller than what is expected 

in the B decay. The term of order A 4 in V cs is necessary to explain K° CP violation. 
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achievement was accomplished using an asymmetric e + e _ collider on the T(4S') which 
was first suggested by Pierre Oddone. The measurements are listed in Table along with 
other previous indications (Groom 2001). 


Table 9. Measurements of sin 2/?. 


Experiment 

sin 2f3 

BABAR 

0.59 ±0.14 ±0.05 

BELLE 

0.99 ±0.14 ±0.06 

Average 

0.79±0.11 

CDF 

o.79i 0 4 L 

ALEPH 

0.84i?;^±0.16 

OPAL 

3.2i^ ± 0.5 


The average value of 0.79 ± 0.11 is taken from BABAR and BELLE only. These two 
measurements do differ by a sizeable amount, but the confidence level that they correctly 
represent the same value is 6%. This value is consistent with what is expected from the 
other known constraints on p and rj. We have 


V 


1 ± 

(1-P) — 


yjl - sin 2 2/3 
sin 2(5 


(99) 


There is a four fold ambiguity in the translation between sin 2/3 and the linear constraints 
in the p — i) plane. These occur at (3, tt/2 — (5, tt + /3 and 3n/2 — (3. Two of these constraints 
are shown in Figure [|T| The other two can be viewed by extending these to negative rj. 
We think rj > 0 based only on measurement of e' in the neutral kaon system. 


7.5 Remarks on Global Fits to CKM parameters 

The fits shown in this paper (Hocker 2001) for p and rj have been done by others with a 
somewhat different statistical framework (Ciuchini 2001) (Mele 1999). The latter group 
uses “Bayesian” statistics which means that they use apriori knowledge of the probability 
distribution functions. The former are termed “frequentists” (Groom 2001), almost for 
the lack of a better term. The frequentists are more conservative than the Bayesians. 

The crux of the issue is how to treat theoretically predicted parameters that are used to 
translate measurements into quantities such as V u b or €k that form constraints in the p — rj 
plane. This of course is a problem because it is difficult to estimate the uncertainties in the 
theoretical predictions. Both groups treat the experimental measurements as Gaussian 
distributions with the la width derived from both the statistical and systematic errors 
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Figure 51. Constraints from sin 2/3 measurement overlaid with other constraints (Hooker 
2001). The inner band is at la while the outer band, shown on one band only, is at 2a. 


combined. Note, that the systematic errors are also difficult sometimes to quantify and 
are not necessarily Gaussian, but they judged to be sufficiently well known as to not cause 
a problem. 

Hocker et al. (Hocker 2001) have decided to use a method which restricts the the¬ 
oretical quantities to a 95% confidence interval where the parameter in question is just 
as likely to lie anywhere in the interval. They call this the “Rfit” method. Of course 
assigning the 95% confidence interval is a matter of judgment which they fully admit. 
Ciuchini et al. (Ciuchini 2001) treat the theoretical errors in the same manner as the 
experimental errors. They call theirs “the standard method” with just a bit of hubris. 
They argue that QCD is mature enough to trust its predictions, that they know the sign 
and rough magnitude of corrections and they can assign reasonable errors, so it would be 
wrong to throw away information. 

Hocker et al. point out an extreme interesting but generally unknown danger with 
the Bayesian approach, which is that in multi-dimension problems the Bayesian treatment 
unfairly predicts a narrowing of possible results. The following discussion will demonstrate 
this. 

Let Xi denote N theoretical parameters over the identical ranges [—A,+A]; then the 
theoretical prediction is 

N 

T { p N) = lixi . ( 100 ) 

i 

In the 95% scan scheme [Tp^] = [—A^+A^] while in the Bayesian approach the con- 
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voluted Probability Density Function (PDF) is 

/ +oo r+oa N , AT \ 

... I] dxiGixiMT-T™) , (101) 

-oo J — oo ■ 

where the G(xi) are PDF’s for each individual variable taken to be equal here. This 
function has a singularity in p(T ) that goes as (— InT) 1 * -1 , so it increases as N increases. 

Now suppose G{xi) is flat, then for N = 1 both approaches are the same, but for 
N > 2, the Bayesian approach gets a p(T) that peaks at zero. In effect, when the number 
of theoretical predictions entering the computation increases, and hence our knowledge of 
the corresponding observable decrease the Bayesian approach claims the converse. 

Lets look at a specific example: here N — 3 and A = \/ 3 . Consider both the sum 
Ts = x\ + X2 + £3 and product Tp = X1X2X3 distributions. For Rfit the allowed ranges are 
identical being [—3A,3A], The left side of Figure |52| shows the probability density p(T) 
for Tg, while the right side shows p(T) for Tp with G(x) in the Bayesian case being either 
a Gaussian with a = 1 (solid lines) or a uniform distribution over the range [— a / 3 , +\/ 3 ] 
(dashed lines). The later distribution is closest to the Rfit method where the allowed 
range for either T is [—3\/3, +3\/3] indicated by the arrows. 

In the Rfit scheme the two predictions for Tp and Tg are identical, while in the Bayesian 
scheme there is large difference in the PDF’s and it really doesn’t matter if a Gaussian or 
uniform distribution is chosen. There is a clear distinction between the Rfit and Bayesian 
predictions for Tp in this case, and the Bayesian one is unreasonable because it produces 
a very narrow PDF peaked very close to zero. 

An example of how this can effect the results is shown on Figure [53] where predictions 
of sin 2/3 from the indirect measurements are shown for the Rfit technique and either 
uniform or Gaussian Bayesian PDF’s. The predictions are quite different. 


8 New Physics Studies 

8.1 Introduction 

There are many reasons why we believe that the Standard Model is incomplete and there 
must be physics beyond. One is the plethora of “fundamental parameters,” for example 
quark masses, mixing angles, etc... The Standard Model cannot explain the smallness of 
the weak scale compared to the GUT or Planck scales; this is often called “the hierarchy 
problem.” In the Standard Model it is believed that the CKM source of CP violation 
extensively discussed here is not large enough to explain the baryon asymmetry of the 
Universe (Gavela 1993); we can also take the view that we will discover additional large 
unexpected effects in b and/or c decays. Finally, gravity is not incorporated. John Ellis 
said “My personal interest in CP violation is driven by the search for physics beyond the 
Standard Model” (Ellis 2000). 
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Figure 52. Convolution of the sum Ts = x\ -\-X 2 + £3 for both Rfit and Bayesian methods 
(left) and the product Tp = X 1 X 2 X 3 (right) of 3 parameters for the Bayesian method only. 
Plotted is the PDF p(T ) obtained using for G(x) a uniform (solid lines, A = \ff>) or a 
Gaussian (dashed lines, a = 1) distribution. Both PDF's ofT P have a singularity at the 
origin which is not shown. The Rfit ranges of Ts and T P are indicated by the arrows 
located in both instances at ±3\/3. From (Flocker 2001). 

We must realize that all our current measurements are a combination of Standard 
Model and New Physics; any proposed models must satisfy current constraints. Since the 
Standard Model tree level diagrams are probably large, lets consider them a background 
to New Physics. Therefore loop diagrams and CP violation are the best places to see New 
Physics. 

The most important current constraints on New Physics models are 

• The neutron electric dipole moment, d n < 6.3 x 10 _26 e-cm. 

. B(b -> sy) = (3.23 ± 0.42) x 10“ 4 and B(b -> sl + t~) < 4.2 x 10~ 5 . 

• CP violation in K L decay, Ck = (2.271 ± 0.017) x 10~ 3 . 

• B° mixing parameter A m c i = (0.487 ± 0.014) ps” 1 . 

8.2 Generic Tests for New Physics 

We can look for New Physics either in the context of specific models or more generically, 
for deviations from the Standard Model expectation. 
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Figure 53. Comparison between Rfit (broad solid curve) and Bayesian fits for the indirect 
CKM constraints on sin 2/3. For the Bayesian fits: Gaussian (uniform) theoretical PDF’s 
are depicted as dashed-dotted (dotted) curves. (This example has been outdated by newer 
data.) 


One example is to examine the rare decays B —> K£ + £~ and B —> K*£ + £~ for branch¬ 
ing ratios and polarizations. According to Grenb et ah, “Especially the decay into K* 
yields a wealth of new information on the form of the new interactions since the Dalitz 
plot is sensitive to subtle interference effects” (Greub 1995). 

Another important tactic is to test for inconsistencies in Standard Model predictions 
independent of specific non-standard models. 

The unitarity of the CKM matrix allows us to construct six relationships. These may 
be thought of as triangles in the complex plane shown in Figure [uf]. (The bd triangle is 
the one depicted in Figure [48].) 



Figure 54. The six CKM triangles. The bold labels, e.g. ds refer to the rows or columns 
used in the unitarity relationship. The angles defined in equation ( \1 OQ ) are also shown. 


All six of these triangles can be constructed knowing four and only four independent 
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angles (Silva 1997) (Aleksan 1994). These are defined as: 


/3 = arg 


X = arg 


v*v« 

v cb v; d/ 

V^Vcb ' 

V t * s v tb/ 


7 = arg 


X = arg 


v: b VuA 

v; b v cd ) 

v: d v U s ' 

V cd V cs , 


( 102 ) 


(103) 


(a can be used instead of 7 or f3.) Two of the phases /3 and 7 are probably large while x 
is estimated to be small ~ 0 . 02 , but measurable, while y' is likely to be much smaller. 

ft has been pointed out by Silva and Wolfenstein (Silva 1997) that measuring only 
angles may not be sufficient to detect new physics. For example, suppose there is new 
physics that arises in B° — B mixing. Let us assign a phase 9 to this new physics. If we 
then measure CP violation in B° —> J/ ip K s and eliminate any Penguin pollution problems 
in using B° —> 7r + 7r~, then we actually measure 2/3 1 = 2/3 + 9 and 2a' = 2a — 9. So while 
there is new physics, we miss it, because 2/3' + 2a' = 2a + 2(3 and a' + (3' + 7 = 180°. 

8.2.1 A Critical Check Using x 


The angle y, defined in equation |102| , can be extracted by measuring the time dependent 
CP violating asymmetry in the reaction B s —> J/iprf'\ or if one’s detector is incapable of 
quality photon detection, the J/ip(p final state can be used. However, in this case there 
are two vector particles in the final state, making this a state of mixed CP, requiring a 
time-dependent angular analysis to extract y, that requires large statistics. 

Measurements of the magnitudes of CKM matrix elements all come with theoretical 
errors. Some of these are hard to estimate. The best measured magnitude is that of 
A = \V us /V ud \ = 0.2205 ± 0.0018. 

Silva and Wolfenstein (Silva 1997) (Aleksan 1994) show that the Standard Model can 
be checked in a profound manner by seeing if: 


sm y = 


K. 


K 


ud 


sin/3 sin 7 
sin(/3 + 7 ) 


(104) 


Here the precision of the check will be limited initially by the measurement of siny, not 
of A. This check can reveal new physics, even if other measurements have not shown any 
anomalies. 


Other relationships to check include: 


sm y = 


V ub 


V cb 

v td 


V, 


ts 


siny sin(/3 + 7 ) 
sin /3 

sin j3 sin(/3 + 7 ) 


(105) 


sm x = 


sm 7 


(106) 
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The astute reader will have noticed that these two equations lead to the non-trivial 
relationship: 


sin 2 /3 


V td 

2 

• 2 

v uh 

Vts 

= sm 7 

v cb 


(107) 


This constrains these two magnitudes in terms of two of the angles. Note, that it is in 
principle possible to determine the magnitudes of \V ub /V cb \ and \V t d/V ts \ without model 
dependent errors by measuring (3, 7 and x accurately. Alternatively, f3, 7 and A can be 
used to give a much more precise value than is possible at present with direct methods. 
For example, once /3 and 7 are known 


Vub 2 _ 2 sin 2 0 
V cb sin 2 (/3 + 7 ) 


(108) 


Table | 8 . 2 . 1 | lists the most important physics quantities and the decay modes that 
can be used to measure them. The necessary detector capabilities include the ability to 
collect purely hadronic final states labeled here as “Hadron Trigger,” the ability to identify 
charged hadrons labeled as “AV sep,” the ability to detect photons with good efficiency 
and resolution and excellent time resolution required to analyze rapid B s oscillations. 
Measurements of cos(20) can eliminate 2 of the 4 ambiguities in 0 that are present when 
sin( 20 ) is measured. 


Physics 

Quantity 

Decay Mode 

Hadron 

Trigger 

AV 

sep 

7 

det 

Decay 

time cr 

sin( 2 a) 

B° — > p7T — ■> 7T + 7r - 7r° 

V 

V 

V 


cos( 2 ci) 

B° — > p-rr — > 7r + 7r _ 7r° 

V 

V 

V 


sign(sin( 2 o:)) 

B° — > pTT & B° — > 7T + 7r~ 

V 

V 

V 


sin( 7 ) 

B s -7 DfK* 

V 

V 


V 

sin( 7 ) 

B~ -»■ D°K~ 

V 

V 



sin( 7 ) 

B° -> 7T+7T- & B s -> K+K~ 

V 

V 


V 

sin( 2 y) 

B s — > J/W, J/'if'q 



V 


sin( 20 ) 

B ° -► J / if K s 





cos( 20 ) 

B° -> J/ifK°, I<° -> 7 t£u 


V 



cos (2 0 ) 

B° -»■ J/ifK*° & A, -> J/00 





x s 

B s ^ D+tt- 

V 



V 

AT for B s 

A, -»■ J/ijrrf, D+n~, K + K~ 

V 

V 

V 

\/ 


Table 10. Required CKM Measurements for b’s 
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8.2.2 Finding Inconsistencies 

Another interesting way of viewing the physics was given by Peskin (Peskin 2000). Non- 
Standard Model physics would show np as discrepancies among the values of (p, rj) derived 
from independent determinations using CKM magnitudes (| V u b/V c b\ and \Vtd/V ts \), or B% 
mixing ((3 and a), or B s mixing (y and 7 ). 

8.2.3 Required Measurements Involving (3 

Besides a more precise measurement of sin 2/3 we need to resolve the ambiguities. There 
are two suggestions on how this may be accomplished. Kayser (Kayser 1997) shows that 
time dependent measurements of the final state J/i^K 0 , where K° —> 7 rtV, give a direct 
measurement of cos(2 (3) and can also be used for CPT tests. Another suggestion is to use 
the final state J/xl>K*°, K*° —> K s n°, and to compare with B s —> Jto extract the 
sign of the strong interaction phase shift assuming SU(3) symmetry, and thus determine 
cos(2/3) (Dighe 1998). 

8.2.4 Required Measurements Involving a and 7 

It is well known that sin(2/3) can be measured without problems caused by Penguin 
processes using the reaction B° —> J /K s . The simplest reaction that can be used to 
measure sin(2o;) is B° —> 7i + 7i~. This reaction can proceed via both the Tree and Penguin 
diagrams shown in Figure [57]. 

Current measurements shown in Table [|show a large Penguin component. The ratio 
of Penguin amplitude to Tree amplitude in the 7 r + 7 r - channel is about 15% in magnitude. 
Thus the effect of the Penguin must be determined in order to extract a. The only model 
independent way of doing this was suggested by Gronau and London, but requires the 
measurement of B T —» and B° —> 7 i°7i°, the latter being rather daunting. 

There is however, a theoretically clean method to determine a. The interference 
between Tree and Penguin diagrams can be exploited by measuring the time dependent 
CP violating effects in the decays B° —> pir —> 7 r + 7 r _ 7 r° as shown by Snyder and Quinn 
(Snyder 1993). 

The pn final state has many advantages. First of all, it has been seen with a relatively 
large rate. The branching ratio for the p°rc + final state as measured by CLEO is (1.5±0.5± 
0.4) x 10 -5 , and the rate for the neutral B final state p ±, K T is (3.51^0 ± 0.5) x 10~ 5 , while 
the p°tt° final state is limited at 90% confidence level to < 5.1 x 10 -6 (Gao 1999). (BABAR 
(Bona 2001) measures B (B° —> p ± tt :¥ ) as (28.9 ± 5.4 ± 4.3) x 10“ 6 .) These measurements 
are consistent with some theoretical expectations (Ah 1999). Furthermore, the associated 
vector-pseudoscalar Penguin decay modes have conquerable or smaller branching ratios. 
Secondly, since the p is spin-1, the 7 r spin-0 and the initial B also spinless, the p is fully 
polarized in the ( 1 , 0 ) configuration, so it decays as cos 2 0 , where 6 is the angle of one of 
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the p decay products with the other 7 r in the p rest frame. This causes the periphery of 
the Dalitz plot to be heavily populated, especially the corners. A sample Dalitz plot is 
shown in Figure [55| This kind of distribution is good for maximizing the interferences, 
which helps minimize the error. Furthermore, little information is lost by excluding the 
Dalitz plot interior, a good way to reduce backgrounds. 



Figure 55. The Dalitz plot for B° —► pn —► n + 7i i r° from Snyder and Quinn. 

To estimate the required number of events Snyder and Quinn preformed an idealized 
analysis that showed that a background-free, flavor-tagged sample of 1000 to 2000 events 
was sufficient. The 1000 event sample usually yields good results for a, but sometimes does 
not resolve the ambiguity. With the 2000 event sample, however, they always succeeded. 

This technique not only finds sin(2a), it also determines cos(2a), thereby removing 
two of the remaining ambiguities. The final ambiguity can be removed using the CP 
asymmetry in B° —> 7 r + 7 r - and a theoretical assumption (Grossman 1997). 

Several model dependent methods using the light two-body pseudoscalar decay rates 
have been suggested for measuring 7 The basic idea in all these methods can be sum¬ 
marized as follows: B° —> 7 r + 7 r - has the weak decay phase 7 . In order to reproduce 
the observed suppression of the decay rate for 7 r + 7 r~ relative to K± 7 t t we require a large 
negative interference between the Tree and Penguin amplitudes. This puts 7 in the range 
of 90°. There is a great deal of theoretical work required to understand rescattering, 
form-factors etc... We are left with several ways of obtaining model dependent limits, 
due to Fleischer and Mannel (Fleischer 1998), Neubert and Rosner (Neubert 1998), Fleis¬ 
cher and Buras (Fleischer 2000), and Beneke et al. . (Beneke 2001). The latter make a 
sophisticated model of QCD factorization and apply corrections. Figure |56| shows values 
of 7 that can be found in their framework, once better data are obtainable. 

I 11 fact, it may be easier to measure 7 than a in a model independent manner. There 
have been two methods suggested. 

(1) Time dependent flavor tagged analysis of B s —> DfK T . This is a direct model 
independent measurement (Du 1986) (Aleksan 1992) (Aleksan 1995). Here the Cabibbo 
suppressed V u & decay interferes with a somewhat less suppressed V c b decay via B s mixing 
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Figure 56. Model predictions from (Beneke 2001) as a function of the indicated 
rate ratios. The dotted curve shows the predictions froin naive factorization. The curved 
bands show the total model uncertainties where the inner band comes from theoretical 
input uncertainties, while the outer band allows for errors to corrections on the theory. 
The specific sensitivity to \V u b\ is showed as the dashed curves. The gray bands show the 
current data with a la error while the lighter bands are at 2a. 


as illustrated in Figure [57| (left). Even though we are not dealing with CP eigenstates 
here there are no hadronic uncertainties, though there are ambiguities. 

(2) Measure the rate differences between B~ —> D°K~ and B + —> D°K + in two 
different D° decay modes such as K~n + and Ii + K~. This method makes use of the 
interference between the tree and doubly-Cabibbo suppressed decays of the D°, and does 
not depend on any theoretical modeling (Atwood 1997) (Gronau 1991). See Figure [57] 
(right). 


8.3 New Physics Tests in Specific Models 

8.3.1 Supersymmetry 

Supersymmetry is a kind of super-model. The basic idea is that for every fundamental 
fermion there is a companion boson and for every boson there is a companion fermion. 
There are many different implementations of couplings in this framework (Masiero 2000). 
In the most general case we pick up 80 new constants and 43 new phases. This is clearly 
too many to handle so we can try to see things in terms of simpler implementations. In the 
minimum model (MSSM) we have only two new fundamental phases. One 9d would arise 
in B° mixing and the other, 6 a , would appear in B° decay. A combination would generate 
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Figure 57. (left) The two diagram diagrams for B s —> DfK T that interfere via B s 
mixing, (right) The two interfering decay diagrams for B~ —* D° K~ where one is ah —> u 
transition and the other a doubly-Cabibbo suppressed decay. 


CP violation in D° mixing, call it (pKn when the D° —> K n + (Nir 1999). Table |lT| shows 
the CP asymmetry in three different processes in the Standard Model and the MSSM. 

Table 11. CP Violating Asymmetries in the Standard Model and the MSSM. 


Process 

Standard Model 

New Physics 

B° - 

+ J/*fK s 

sin 2/3 

sin 2 {(3 + 0 D ) 

B° - 


sin 2(3 

sin 2 {(3 + 0 D + 6 A ) 

D° - 

-> K~tt+ 

0 

~ sin f) Kn 


Two direct effects of New Physics are clear here. First of all, the difference in CP 
asymmetries between B° —> .J / ipK s and B° —> <f>K s would show the phase 4>a- Secondly, 
there would be finite CP violation in D° —> K~n + where none is expected in the Standard 
Model. 


Manifestations of specific SUSY models lead to different patterns. Table |12| shows the 
expectations for some of these models in terms of these variables and the neutron electric 
dipole moment d^\ see (Nir 1999) for details. Note, that “Approximate CP” has already 
been ruled out by the measurements of sin 2/3. 


In the context of the MSSM there will be significant contributions to B s mixing, 
and the CP asymmetry in the charged decay B T —> cf)K T . The contribution to B s 
mixing significantly enhances the CP violating asymmetry in modes such as B s —> J / ifrj. 
(Recall the CP asymmetry in this mode is proportional to sin2y in the Standard Model.) 
The Standard Model and MSSM diagrams are shown in Figure The expected CP 
asymmetry in the MSSM is ~ sin^coS(^ 4 Sin(Am s i), which is approximately 10 times 
the expected value in the Standard Model (Hinchliff 2001a). 
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Table 12. Some SUSY Predictions. 


Model 

d N x 10 25 

0 D 

0 A 

sin (pKn 

Standard Model 

< 10” 6 

0 

0 

0 

Approx. Universality 

> 10~ 2 

0(0.2) 

0(1) 

0 

Alignment 

> 10~ 3 

C>(0.2) 

0(1) 

0(1) 

Heavy squarks 

~ 10" 1 

0(1) 

0(1) 

O(10~ 2 ) 

Approx. CP 

~ 10” 1 

-P 

0 

O(10~ 3 ) 



*** _L _L 

w + h 2 + 


Figure 58. The Standard Model (left) and MSSM (right) contributions to B° mixing. 


We observed that a difference between CP asymmetries in B° —> ,J/ ip K s and <f>K s arises 
in the MSSM due to a CP asymmetry in the decay phase. It is possible to observe this 
directly by looking for a CP asymmetry in B T —> <f>K T . The Standard Model and MSSM 
diagrams are shown in Figure 59. Here the interference of the two diagrams provides the 
CP asymmetry. The predicted asymmetry is equal to (Mw / msquark) 2 sin^ in the MSSM, 
where m squar k is the relevant squark mass (Hinchliff 2001a). 


W _ x 



Figure 59. The Standard Model (left) and MSSM (right) contributions to B —> <fK . 

The <pK and <f>K* final states have been observed, first by CLEO (Briere 2001) and 
subsequently by BABAR (Aubert 2001). The BABAR data is shown in Figure |Hf]. The 
average branching ratio is B(B~ —> <fK~) = (6.8 ± 1.3) x 10~ 6 showing that in principle 
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large samples can be acquired especially at hadronic machines. 



5.20 5.25 5.30 5.20 5.25 5.30 

m ES (GeV/c 2 ) 

Figure 60. Projections of the maximum likelihood fit onto the B candidate mass rn^s 
from BABAR. (a) B + -> <f)K+; (b) B° -> <f>K°; (c) B + -> <t>K*+; (d) B° -> <fK*°. In 
(c) the histogram is the sum of the two contributing K* channels while the shaded area is 
K°tt + alone (the other channel is K + tt°). The solid line shows the projection of the signal 
plus background fit while the dashed line shows the projection of the background only. 


8.3.2 Other New Physics Models 

There are many other specific models that predict New Physics in b decays. I list here a 
few of these with a woefully incomplete list of references, to give a flavor of what these 
models predict. 

• Two Higgs and Multi-Higgs Doublet Models- They predict large effects in €k and CP 
violation in D° —> K~n + with only a few percent effect in B° (Nir 1999). Expect 
to see 1-10% CP violating effects in b —> s-y (Wolfenstein 1994). 

• Left-Right Symmetric Model- Contributions compete with or even dominate over 
Standard Model contributions to B^ and B s mxing. This means that CP asymme¬ 
tries into CP eigenstates could be substantially different from the Standard Model 
prediction (Nir 1999). 

• Extra Down Singlet Quarks- Dramatic deviations from Standard Model predictions 
for CP asymmetries in b decays are not unlikely (Nir 1999). 
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• FCNC Couplings of the Z boson- Both the sign and magnitude of the decay leptons 
in B —> K*£ + £~ carry sensitive information on new physics. Potential effects are 
on the of 10% compared to an entirely negligable Standard Model asymmetry of 
~ 10~ 3 (Buchalla 2000). These models also predict a factor of 20 enhancement of 
b —> d£ + £~ and could explain a low value of sin 2/3 (Barenboim 2001a). 


• Noncommutative Geometry- If the geometry of space time is noncommutative, i.e. 
[Xfj,,x v ] = i0 llu , then CP violating effects may be manifest a low energy. For a scale 
<2 TeV there are comparable effects to the Standard Model (Hinchliffe 2001b). 


• MSSM without new flavor structure- Can lead to CP violation in b —> sy of up to 5% 
(Bartl 2001). Ali and London propose (Ah 1999) that the Standard Model formulas 
are modified by Supersymmetry as 


A m d 
A m s 

I ex | 


Am d (SM) 1 + / (m x ±, m iR , m H ±, tanfty 
Am s (SM) 1 + / (jn x ±, m ^ R , rrifj ±, tanfi^j 

GpfflMKM^y 2 ^ 6 —\ \ t t ( \ \ 

6V2tt 2 A M k Bk ^ A A ^ ct ^ yciyt ’ ~ Vcc ) 
+Vttytfs(yt) [l + / (m x ±, m iR , m H ±,tan{. 3)] A 2 A 4 (1 - p) 


(109) 

( 110 ) 


( 111 ) 


where A m(SM) refers to the Standard Model formula and the expression for \ek\ 
would be the Standard Model expression if / were set equal to zero. Ali and London 
show that it is reasonable to expect that 0.8 > / > 0.2 so since the CP violating 
angles will not change from the Standard Model, determining the value of (p, rf) 
using the magnitudes Am s /Amd and \ck\ will show an inconsistency with values 
obtained using other magnitudes and angles. 


• Extra Dimensions- We are beginning to see papers predicting b decay phenom¬ 
ena when the world has extra dimensions. See (Agashe 2001) (Barenboim 2001b) 
(Branco 2001) (Chang 2001) and (Papavassiliou 2000). 


I close this section with a quote from Masiero and Vives (Masiero 2001): “The rele¬ 
vance of SUSY searches in rare processes is not confined to the usually quoted possibility 
that indirect searches can arrive ‘first’ in signaling the presence of SUSY. Even after the 
possible direct observation of SUSY particles, the importance of FCNC and CP violation 
in testing SUSY remains of utmost relevance. They are and will be complementary to 
the Tevatron and LHC establishing low energy supersymmetry as the response to the 
clectroweak breaking puzzle.” 

I agree, except that 1 would replace “SUSY” with “New Physics.” 
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8.4 Expected Data Samples 

It is clear that precision studies of b decays can bring a wealth of information to bear on 
new physics, that probably will be crucial in sorting out anything seen at the LHC. This 
is possible because we do expect to have data samples large enough to test these ideas 
from existing and approved experiments. In Table [L3| I show the expected rates in BTeV 
for one year of running (10 7 s) and an e + e~ H-factory operating at the T(4S') with a total 
accumulated sample of 500 fb—1, about what is expected around 2006. (LHCb numbers 
are the same order of magnitude as the BTeV numbers for many of the modes.) 


Table 13. Comparisons of BTeV and B-factory Yields on Different Time Scales. 


Mode 

BTeV (10 7 )i 
Yield Tagged' 1 ' 

S/B 

R-factory (500 fb *) 
Yield Tagged^ S/B 

B s - 


22000 

2200 

>15 




B~ -> 

4>K~ 

11000 

11000 

>10 

700 

700 

4 

B° -> 

f>K s 

2000 

200 

5.2 

250 

75 

4 

B° -»• 

K*°n+ir 

4400 

4400 

11 

~50 

~50 

? 

D*+ - 

-> 7r + P°; D° —> K~ir + 

~ 10 8 

~ 10 8 

large 

8 x 10 5 

8 x 10 5 

large 


f Tagged here means that the intial flavor of the B is determined. 


9 Conclusions 

I have attempted to cover the length and breadth of b physics, and have only scratched 
the surface. There is much more to be said and much more to learn. Why are there three 
families? What is the connection with neutrinos and that mixing matrix? How do we 
explain the mystery of flavor? These and many more unanswered questions I leave to the 
students. 
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